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In the first series of measurements the source was one 
No. 6 blasting cap, while in the second series one pound of 
50 percent nitroglycerine stick dynamite was used as 
source. A telephone carbon button microphone was the 
receiver and was held fixed in location while shots were 
fired successively at 5, 8, 12.5, 25, 35, 50, 100 and 600 
meter distances. A two-element string oscillograph was 
used for timing, one element recording the instant of firing, 
the other element recording the arrival of the wave at the 
microphone. The work was done at time of no perceptible 
wind; air temperatures were measured carefully; no humid- 
ity measurements were made. Travel times could be read 


reliably to 10-* second and distance measurements were at 
least correspondingly good. Instantaneous speeds were ob- 
tained by plotting time computed minus time observed 
against distance and measuring slopes of the resulting 
curve. Since this work was incidental to seismic prospect- 
ing, the observations were not quite as numerous as those 
of von Angerer and Ladenburg. However, the results are 
similar as regards abnormally high speeds near the source. 
Also, the use of instantaneous speeds appears to show abnor- 
mall ylow speeds alittle farther from the source, perhaps masked 
in the work of von Angerer and Ladenburg due to the use 
by them of average speeds instead of instantaneous speeds. 





I. INTRODUCTION 


NE of us, from experience with seismic 

prospecting, and the other,! from experience 
with echo sounding and sound ranging in water, 
had long held the opinion that sound waves of 
finite amplitude are propagated with abnormally 
high speeds. Since for our experimental work in 
refraction seismic prospecting we have at hand 
adequate means of measuring travel times, it 
was decided to carry out two series of measure- 
ments of travel times in air over short distances, 
one series with one blasting cap as source, the 
other series with one pound of 50 percent nitro- 
glycerine stick dynamite as source. The results 
of the measurements, which were carried out in 
October, 1933, are as stated in this paper. Not 
much search of the literature for accounts of 
similar work was made prior to the carrying out 
of the observations and the computation of the 
results. Nor has an exhaustive search yet been 
made. However, the work of von Angerer and 


1 Jerry H. Service, The Transmission of Sound through 
Sea Water, J. Frank. Inst. 206, 802 (1928). 


Ladenburg? is the only published work that has 
come to our knowledge that was carried out 
under conditions comparable to those involved 
in this work of ours. Their observations at short 
distances were a little more numerous than ours. 
Our time measurements, because of a larger 
time scale on the oscillograph records, may have 
been a little more accurate than theirs. They 
computed average speeds over various sections 
of the course ; we computed instantaneous speeds 
at various distances from the source. Their 
results showed abnormally high speeds near the 
source, decreasing asymptotically to the normal 
value (for faint sounds) within a comparatively 
short distance. Our results are consistent with 
theirs, except that we found abnormally high 
speeds near the source, abnormally low speeds a 
little farther from the source, then normal speed 
for greater distances. We cannot but wonder 
whether von Angerer and Ladenburg would not 
find these abnormally low speeds also, if they 
would compute instantaneous, rather than aver- 
age, speeds from their data. 

2E. von Angerer and R. Ladenburg, Experimentelle 


Beitraege zur Ausbreitung des Schalles in der freien Atmos- 
phaere, Ann. d. Physik 66, 293 (1921). 
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TABLE I. Data using one No. 6 blasting cap. 
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TABLE II. Data using 1 1b. 50 percent nitroglycerine. 











x vo x v0 
No. (meters) (°C (m/sec.) <x/vo T AT Est No. (meters) (°C (m/sec.) x/v0 rs AT Est. 
1 5.005 11.6 337.9 0.0148 0.0079 —0.0069 +2 g& 5.005 11.1 337.6 0.0148 0.0078 —0.0070 +3 
2 8.007 11.6 337.9 0.0237 0.0165 —0.0072 4 7 8.007 11.0 337.5 0.0237 0.0146 —0.0091 2 
3 12.506 11.6 337.9 0.0370 0.0242 —0.0128 2 6 12.506 11.3 337.7 0.0370 0.0257 —0.0113 2 
4 25.000 11.6 337.9 0.0740 0.0632 —0.0108 2 5 25.000 12.1 338.2 0.0739 0.0596 —0.0143 2 
9 35.305 9.9 336.8 0.1048 0.0928 —0.0120 4 12 35.305 10.3 337.1 0.1047 0.0887 —0.0160 2 
10 50.003 10.8 337.4 0.1482 0.1402 —0.0080 3 11 50.003 10.4 337.1 0.1483 0.1307 —0.0176 3 
14 100.00 9.5 336.6 0.2971 0.2914 —0.0057 2 13 66.187 8.9 336.2 0.1968 0.1838 —0.0130 2 
16 158.44 8.0 335.6 0.4721 0.4705 —0.0016 6 15 100.00 10.0 336.9 0.2968 0.2836 —0.0132 3 
17 600.03 6.0 334.4 1.7943 1.7860 —0.0083 10 18 600.03 5.8 334.3 1.7949 1.7800 —0.0149 6 








II. EXPERIMENTAL PROCEDURE 


An ordinary telephone transmitter carbon 
button microphone was used as receiver. It was 
supported by means of rubber bands at the 
center of a hoop about 18 inches in diameter so 
that it faced the line of shot points and was 
about four feet above the ground. The micro- 
phone was connected through a 3-stage vacuum 
tube amplifier to one element of a General Radio 
string oscillograph. 

Stakes about four feet high were driven into 
the ground along a straight line through the 
microphone, at distances approximately 5, 8, 
12.5, 25, 35, 50, 100 and 600 meters from the 
microphone. The distances were then taped 
with a 100-foot steel tape. The charges were 
suspended loosely from the tops of the stakes. 

The microphone was kept in a fixed location 
and the charges were fired successively at the 
various distances. In the first series the charge 
was a No. 6 electric blasting cap fired at each 
distance; in the second series the charge was 


two sticks (approximately one pound total) of 
50 percent nitroglycerine powder fired at each 
distance by means of an electric blasting cap. 
Instant of firing was recorded on the second 
element of the string oscillograph by wrapping 
a wire around the blasting cap, connecting this 
wire in series with a dry cell and the primary of 
a transformer, and connecting the secondary of 
the transformer to the galvanometer string 
through a 3-stage amplifier alike in all respects 
to the one used for the microphone. The charges 


- were fired at times when no motion of the air 


sufficient to move an anemometer was to be 
noted. Air temperatures near the ground were 
measured at each shot by means of standardized 
thermometers. No humidity measurements were 
made. No barometer readings were made, but 
barograph records are available at the college 
for the time of the shots. 

Fig. 1 shows some of the records. Travel times 
were scaled off independently by the two of us 
and averaged, with the results given in Tables I 
and IT. Our agreement with each other and with 





(b) 


Fic. 1. (a) Record No. 6. 1 lb. dynamite, 125 meters. Note arrival of echo at 0.151. (b) Record No. 3. 1 blasting cap, 
125 meters. 
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Fic. 2. Experimental relation between AT (travel time 
observed minus travel time computed) and x (distance 
from disturbance to receiver) Curve A for No. 6 blasting 
cap and Curve B for one pound 50 percent nitroglycerine 
powder. 


ourselves at different times indicates that the 
travel times are probably generally correct to 
10-* second. Taking into consideration slight 
uncertainties on some of the records as to just 
where on the record the reception of the wave 
is represented, however, the last columns of 
Tables I and II give rather pessimistic estimates 
of the uncertainties of the travel times. The 
long time marks on the records are one-fiftieth 
second apart, very nearly, and the short time 
marks are at one-tenth second intervals, very 
nearly. The timing apparatus is controlled by a 
tuning fork and was standardized at several fork 
temperatures against a Bond break circuit 
chronometer, which was checked daily (during 
the standardization) against radio time signals 
from Arlington. The temperature of the tuning 
fork of the time marking apparatus was recorded 
at each shot, as measured by a calibrated 
thermometer. 

Referring to Tables I and II, the first column 
gives the order of firing ; column two, x in meters, 
records the distance from the charge to the 
microphone; é is the air temperature between 
charge and microphone; wv is the accepted® 
normal speed of sound in air at temperature f 
given by the equation 1 =330.6+0.63t; x/vo is 
self-explanatory, being the normal travel time 
for distance x; T is the observed travel time in 
seconds. AT is T minus x/v, while the last 
column of each table gives the uncertainty in 
the fourth decimal place of T and AT as already 
explained. 


*T. Vautier, Acoustique—Sur la propogation d’ondes 
aériennes produites par des étincelles ou des amorces, 
Comptes rendus 181, 1055 (1925). 
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Fic. 3. A portion of Fig. 2 drawn to larger scale, 


In Figs. 2 and 3 are shown AT plotted as 
ordinate against x as abscissa, Fig. 3 being 
merely a large scale reproduction of the part of 
Fig. 2 that involves the shorter distances and 
the abnormal speeds. 


III. CompuTATION OF INSTANTANEOUS SPEEDS 
Defining x, T and wv as above gives, 


T=x/vo+ AT. (1) 


Now the instantaneous speed, v, is given by 
v=dx/dT, and by differentiating (1) we have 
dT /dx=1/v9+d(AT) /dx, 
1 


v= . 3 
1/v9+slope of AT=f(x) curve ( 


Three points, at 5, 12.5 and 158 meters, on 
the A curves of Figs. 2 and 3 were disregarded 
in drawing the curves, because the authors could 
not draw smooth curves through these points. 
No points were disregarded on curves B. 

Slopes of the AT =f(x) curves in Fig. 3 were 
measured by the two of us independently and 
averaged, with results as recorded in the second 
columns of Tables III and IV. From these 
slopes 1/v and v were computed with the results 
given in the third and fourth columns of Tables 
III and IV and shown graphically in Fig. 4. 
The values of v given below Tables III and IV 
are the accepted values for the average air 
temperatures of first and second series of shots, 
respectively. In other words, the average of the 
air temperatures in Table I is 10.1, and likewise 
for Table IT. 


(2) 


whence 
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TABLE III. Computed data (blasting caps). 
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TABLE IV. Computed data (1 lb. 50 percent nitroglycerine). 








Slope of AT =f(x) 


Slope of AT =f(x) 

















x (meters) curve 1/e v (m/sec.) x (meters) curve 1/v v (m/sec.) 
5 —0.658 K 10-3 2.309 K 1074 433.1 5 —0.79 X1073 2.18 10-3 459. 
YP —0.436 XK 1078 2.531 10-3 395.1 7.5 —0.62 X10-3 2.35 1073 426. 
10 —0.337 X 107" 2.630 X 1073 380.2 10 —0.469 X 1073 2.498 K 1073 400.3 
15 —0.245 X 1073 2.722 X10-3 367.4 12.5 —0.383 XK 1073 2.584 x 107% 387.0 
20 —0.176 X 1074 2.791 K 10-3 358.3 15 —0.318 KX 10-3 2.649 X 107% 377.5 
25 —0.144 «107% 2.823 X 1073 354.2 20 —0.230 X 107% 2.737 X 10-3 365.4 
30 —0.091 K 10-4 2.876 X1073 347.7 25 —0.193 K 1073 2.774 X 1073 360.5 
35 +0.005 K 1073 2.972 X 10-3 336.5 30 —0.167 X 1073 2.800 K 10-3 357.1 
40 +0.146 X 1073 3.113 K 10-3 321.2 40 —0.100 X 10-3 2.867 X 107-3 348.8 
50 +0.52 K1073 3.49 «1073 287. 45 —0.064 X 1073 2.903 X 10-3 344.5 
58 —0.005 K 1073 2.962 K 1073 *337.6 50 —0.012 X1073 2.955 X 1073 338.4 

50.8 0.000 « 10-3 2.967 X 1073 337.0 

58 +0.515 K 1073 3.482 K 1073 287.2 
t=10°.1C. 1% =337.0 and 1/2 =0.002967 65 —0.004 x 107% 2.960 X 1073 *337.5 
* 337.6 m/sec. reduced to 0°C =331.2 m/sec. 

t=10°.1C. 1 =337.0. 1/vo =0.002967 


Fig. 4 summarizes the results of the experi- 
ments. It will be seen that with a No. 6 blasting 
cap as source the abnormally high speed of 
433. m/sec. was found at 5 meters from the 
source ; the speed decreased more or less regularly 
to a normal value at a distance of 35 meters 
from the source. From 35 meters to 58 meters 
the speed was abnormally low. At distances 
greater than 58 meters the speed was about 
normal. 

With the heavier charge (two sticks of dyna- 
mite) a higher speed (459 m/sec.) was measured 
at 5 meters from the source, a somewhat greater 
distance (50.8 meters) was required for decrease 
to normal speed, the same minimum abnormally 
low speed (287 m/sec.) was noted as with the 
blasting caps, but at 58 meters instead of 50 
meters. Normal speed was reached finally with 
the heavier charge at 65 meters instead of at 
58 meters as with the caps. 

The authors do not attempt to explain either 
the abnormally high speeds or the abnormally 
low speeds found as the results of these experi- 
ments. Considering the probable accuracy of the 
time, distance and slope measurements, and the 
regularity of the results, it does not seem likely 
that the results are merely accidental. Possibly 
the high speeds near the source are sound speed 
plus explosion gas speed, and possibly the low 
speeds a little farther out are sound speed minus 
explosive gas recoil speed. 

It should be stated again that confusion of the 
air wave with a wave through the earth was 
guarded against by suspending the charges .in 
the air by light string and by suspending the 
microphone in the air on rather light rubber 
bands. Also, those who have done seismic work 


* 337.5 m/sec. reduced to 0°C =331.1 m/sec. 


over comparatively short distances will realize 
that the arrival of the air wave records itself on 
the oscillogram in a rather characteristic manner, 
not at all likely to be confused with the record 
made by an earth wave that must pass through 
a foot of light string and a foot or so of light 
rubber band. If anyone interested should wish to 
study the oscillograms first-hand they are avail- 
able for that purpose. 

The authors regret that they were not able to 
carry out more numerous observations before 
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Fic. 4. Instantaneous speed of sound in air relative to 
the ground with blasting cap (Curve A) and one pound 50 
percent nitroglycerine powder (Curve B) as sources of 
disturbance. 
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publication. Since their major research duty is 
seismic prospecting, it did not seem feasible to 
spend more field time upon the subject. The work 
is published in the hope that it may stimulate 
further work along the same line on the part of 
other experimenters. It should be stated that 
previous to the work for which data are here 
given the authors carried out other similar 
measurements. The work was not included in 
this report because it was of a somewhat lower 
order of accuracy. The absolute values of speed 
agreed only roughly with those here published ; 
however, similar effects with increasing distances 
from the source, abnormally high speeds, then 
speeds somewhat below normal, then normal 
speed, were found as in the case of the work 
here reported. 

It is interesting to note that extrapolating 


back to zero distance, Fig. 4, gives speeds at the 
source of 516 and 544 meters per sec. for cap 
and dynamite respectively. If the curves of 
Figs. 2 and 3 were to continue back smoothly 
from 5 to zero meters, as they must necessarily 
do in order to arrive at the origin, Fig. 2, the 
extrapolation is evidently legitimate and the 
instantaneous speed of the sound relative to the 
earth would not be less than the above figures. 
As a matter of fact, the average speeds for 5.005 
meters are 633.5 and 641.7 meters per sec. for 
cap and dynamite respectively, and at 5.005 
meters the instantaneous speeds are 433.1 and 
459 meters per sec. making the necessary initial 
speed at least 834 and 825 meters per sec., 
respectively. 

The authors hope to be able to carry out 
similar measurements in water in the near future. 
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Shot Effect and Thermal Agitation in an Electron Current Limited by Space Charge 


G. L. Pearson, Bell Telephone Laboratories, New York 
(Received September 28, 1934) 


Experimental evidence is presented which shows that shot effect is decreased by space charge 
and that the internal impedance of a vacuum tube is subject to thermal agitation. 





INTRODUCTION 


HE space current in a thermionic vacuum 

tube is not a steady flow of electricity but is 
subject to minute irregular fluctuations. The 
two most fundamental causes for these fluctua- 
tions are the random distribution of instants of 
emission of the individual electrons and the 
distribution of these electrons in velocity. The 
random emission produces shot noise which may 
be reduced by the space charge surrounding the 
cathode, while the velocity distribution produces 
thermal noise and is dependent upon the tem- 
perature of the cathode. 

Although plausible theories of these effects 
have been given they have never been checked 
by accurate experiments because of the difficul- 
ties involved. By using two electrode tubes 
capable of producing a large space charge such 
measurements have now been made and are re- 
ported in this paper. 


SHot EFFECT 


The theory of the shot effect in the absence of 
space charge has been studied quite completely 
both theoretically and experimentally and more 
recently Llewellyn! has extended the theory to 
cases where partial temperature saturation exists. 
The general equation developed by him is as 


follows: 
oO 2 co 
E3=2j(—) f LG*df. (1) 
oj 0 


E. is here the mean square shot voltage across 
the output measuring device, 7 the space current, 
j the total current emitted by the cathode, e the 
electronic charge, Z the impedance of the parallel 
combination consisting of the internal plate im- 
pedance, r,, of the vacuum tube and the external 
load impedance 2, G the voltage amplification 





'F, B. Llewellyn, Proc. 1.R.E. 18, 243 (1930). 


between 2» and the output measuring device, and 
df the differential of frequency. 

Eq. (1) shows that as long as space charge is too 
small to affect the flow of current, the mean 
square shot voltage is directly proportional to the 
space current. As the emission is increased, how- 
ever, space charge begins to control, and finally 
limits, the space current so that the value of 
0i1/0j approaches zero. At the same time the in- 
ternal plate impedance is reduced and Z is thus 
decreased. The shot voltage therefore increases 
less rapidly as space charge becomes effective and 
finally decreases rapidly toward zero as complete 
space charge control is reached. 

The effect of space charge on shot noise has 
been demonstrated in a qualitative manner by 
several workers.? A precise experimental verifica- 
tion of Eq. (1) is very difficult because most 
vacuum tubes are incapable of producing good 
temperature saturation and even after obtaining 
saturation an exact measurement of 07/0] is al- 
most impossible. This latter difficulty is particu- 
larly objectionable with thoriated and oxide 
coated cathodes. Accordingly, in the measure- 
ments to be described here, a tube having a pure 
tungsten filament was used. Although the tube 
contained three electrodes the grid and plate were 
tied together in order to simplify the system. 
The plate voltage was kept as low as possible in 
order to eliminate extraneous noise due to ioniza- 
tion of the residual gas in the tube. Under these 
operating conditions the thermal noise in the 
plate circuit of the tube is so small that it may be 
neglected. 


2A. W. Hull and N. H. Williams, Phys. Rev. 25, 147 
(1925); J. B. Johnson, Phys. Rev. 26, 71 (1925); N. H. 
Williams and E. W. Thatcher, Phys. Rev. 39, 474 (1932); 
and reference 1 above. E. W. Thatcher, Phys. Rev. 40, 114 
(1932) has reported experimental evidence to show that 
the mean square shot voltage in the presence of space 
charge is proportional to j(d1/0j)?. He did not, however, 
make corrections for Schottky emission in the determi- 
nation of j. 
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Fic. 1. Noise measuring circuit consisting of experimental 
tube, coupling circuit, amplifying system, thermocouple 
and microammeter. 


The experimental arrangement is shown dia- 
grammatically in Fig. 1. The amplifying system 
and the output measuring device are the same as 
those described in a previous paper.* The emis- 
sion in the experimental tube was varied by 
changing the filament heating current while the 
plate voltage remained constant at 9.4 volts. The 
plate battery was varied to compensate for volt- 
age drop in the load resistance, 7, which was 
1000 ohms. The output deflection (proportional 
to the mean square voltage fluctuations across 1) 
was recorded for widely different space currents. 
These are given in Table I and also plotted as ex- 
perimental points in Fig. 2. It is seen that the 
noise varies qualitatively as previously described 
except that instead of reaching the value zero it 
passes through a minimum value, at about 5 


TABLE I. Shot noise data.* 








E2—Arbitrary 








units 
1 jo 
(m.a.) (ohms) (m.a.) 01/9}, Cal. Obs. 
1.0 995 1.0 1.03 27 25 
2.0 989 1.9 1.03 51 50 
2.5 985 2.4 1.03 65 64 
3.0 978 2.9 1.03 77 77 
3.3 969 3.2 1.03 83 82 
3.5 960 3.4 1.03 87 85 
3.7 944 3.6 1.00 83 84 
3.2 935 3.7 0.95 75 80 
3.9 906 . 3.8 0.90 65 68 
4.0 874 4.0 0.85 57 56 
4.1 849 4.2 0.80 50 45 
4.4 760 4.5 0.65 29 32 
4.6 690 4.8 0.47 13 28 
4.8 620 5.5 0.34 6.4 26 
5.0 585 6.2 0.22 2.6 26 
5.2 562 7.4 0.14 +2 29 
5.5 533 10.5 0.068 0.35 40 








* Those values given in italics were interpolated from the measured 
data. 


3G. L. Pearson, Physics 5, 233 (1934). 
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Fic. 2. Shot effect in the presence of space charge. E 
is the mean square noise voltage across the output meas- 
uring device expressed in arbitrary units. The variation in 
space current was obtained by changing the cathode 
temperature, the plate voltage remaining constant. The 
plotted points represent measured noise, the solid curve 
was calculated by means of Eq. (1), and the dashed curve 
calculated according to the relation suggested by Moullin 
and Ellis. 


milliamperes, and then increases. This latter 
effect is well known and is almost certainly a 
fluctuation in the space charge produced by 
positive ions emitted from the tungsten filament.* 
The experiments can thus distinguish between 
Eq. (1) and the pure shot effect equation over the 
range of currents between the values where space 
charge just becomes perceptible and where the 
space charge is so dense that it can be affected 
appreciably by the positive ions. 

For each value of 2, the total emission of the 
cathode was determined in the following manner. 
Keeping the filament temperature constant at 
that value which gave a fixed space current i 
when the plate voltage was 9.4 volts, the voltage 
in the plate circuit was increased and the values 
of space current recorded. The logarithm of the 
space current was plotted against the square root 
of the plate voltage. At the higher voltages where 
the variation in space current results solely from 
the “Schottky current”’ the plot is a straight line.‘ 
By extrapolation of this line to zero voltage the 
total emission 7) at zero field is obtained. Al- 
though this jo is slightly different from the 7 in 
Eq. (1), jo will be used in the calculations since it 
is impossible to determine the field at the cathode 


*K. T. Compton and I. Langmuir, Rev. Mod. Phys. 2, 
123 (1930). 
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under operating conditions. The error introduced 
by this substitution is very small. The space cur- 
rent 7 was now plotted as a function of jp and 
0i1/j) was found by measuring the slope of this 
curve at the various points. The values of j) and 
01/8jo so obtained are given in Table I. The corre- 
sponding values of the internal tube impedance r, 
were also measured in the order to determine Z, 
the impedance of r, in parallel with ro, as given in 
column 2 of Table I. The values of the impedance 
of the tube used in calculating Z were obtained 
simply from the slope of the V-7 curve, V being 
the plate voltage, and may be expressed by 


( oi )+ ai dj ° 
Nav], ajav 

In Llewellyn’s derivation of Eq. (1) it is assumed 
that 





r,=(d1/dV);. 


Using the former value of r, makes Z at the most 
one percent too large. 

By substituting these quantities in Eq. (1) and 
assigning to G the value required to make the ob- 
served and calculated shot noise identical when 
the space current is 3 milliamperes, the values of 
shot noise given in column 5 of Table I were cal- 
culated. These are also shown in the solid curve 
in Fig. 2. The check between theory and experi- 
ment is remarkably good over the range in which 
the comparison is valid. 


THERMAL AGITATION 


Thermal agitation in the internal plate im- 
pedance of a vacuum tube has been investigated 
theoretically by a number of workers.® It has 
been deduced that the resistive component of 
this impedance should produce thermal noise as 
if it were at the temperature of the cathode and 
that the total thermal noise of the circuit con- 
sisting of the tube impedance in parallel with the 
load resistance should be given by the expression : 


ry 
a \(— ) f Gf. (2) 
Tofp/ \rotrp 


E;* is here the mean square thermal agitation 
voltage across the output measuring device, k is 
Boltzmann’s gas constant, and 7) and Ty are, 





E7?=4k Taro( 1+ 


5 For example see reference 1. 


PEARSON 


respectively, the temperature of the load re- 
sistance and of the cathode expressed in degrees 
Kelvin. 

In order to measure this effect experimentally 
the noise from other causes such as shot effect, 
positive ions, etc., must be kept very low. Ac- 
cordingly the tubes in the measurements de- 
scribed here were worked at high cathode tem- 
peratures in order to insure a high degree of space 
charge at all times and the plate potentials were 
kept low to prevent ionization from electron im- 
pact. Tubes having both oxide coated and thori- 
ated tungsten cathodes were used for the meas- 
urements. Pure tungsten cathodes could not be 
used because of the positive ion emission noted 
above. Although the tubes contained three elec- 
trodes the grid and plate were tied together in 
order to simplify the system. 

The experimental arrangement is the same as 
that used in the shot measurements, as shown in 
Fig. 1. With a fixed filament voltage, correspond- 
ing to a high filament temperature, the noise was 
measured as the internal plate impedance was 
varied by changing the plate voltage. The ampli- 
fying system was calibrated by measuring the 
thermal noise as the load resistance alone was 
varied, with the experimental tube in position 
but with zero plate battery so that its impedance 
was infinite. The noise was the same whether the 
filament was hot or cold. This gave the linear 
relationship shown by the dashed line in Fig. 3 
where ordinates represent mean square noise 
voltage across the output measuring device and 
abscissae represent the resistance of the tube and 
the load in parallel. From these measurements 
the value of G was obtained. The load resistance 
was now fixed at 45,500 ohms and different plate 
voltages were applied to the tube. With increased 
voltage the space current increased while the 
tube impedance and the noise both decreased. 
The noise data obtained with two different tubes 
having oxide coated filaments and one tube hav- 
ing a filament of thoriated tungsten are shown in 
Fig. 3. Curves corresponding to the experimental 
points were calculated by Eq. (2) letting To be 
300 degrees and 7; be for the oxide coated fila- 
ments 650 degrees and for the thoriated filament 
1200 degrees Kelvin. 

It is seen that the experimental points fit the 
theoretical curves for thermal noise very well if 
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PLATE CIRCUIT RESISTANCE IN OHMS 


Fic. 3. Thermal noise in the plate circuit of a vacuum 
tube. The dashed curve represents the thermal noise in the 
load resistance alone. With the tube in the circuit the ex- 
perimental points are represented by: + thoriated tungsten; 
open circle, oxide coated No. 1; and closed circle, oxide 
coated No. 2. The solid curves were calculated by means of 
Eq. (2), Ty having the values shown. 


the filament temperatures are assumed to be 
those stated above. In reality it is known that the 
temperature of the coated filament was about 
1100 degrees Kelvin while that of the thoriated 
filament was approximately 1850 degrees. The 
curve corresponding to this latter temperature 
has been calculated and is shown in the upper 
curve of Fig. 3 for comparison. The cause for this 
discrepancy between theory and experiment is 
not understood. Among the various suggestions 
two may be mentioned.® (1) Of the three velocity 
components of the electrons corresponding to 
temperature not all are equally effective in pro- 
ducing thermal noise. (2) In the thermodynamic 
argument! which leads to Eq. (2) the circuit in- 
cluding filament, plate, and external load are 
first supposed to be maintained at the tempera- 
ture of the cathode. Later the external load is 


° The first suggestion was made by Dr. Llewellyn, the 
second by Dr. Johnson. 
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permitted to take any other temperature. In the 
experimental circuit, however, the plate is more 
nearly at the lower temperature of the external 
load than at that of the cathode. Because of the 
close electrical reaction between the plate and the 
space charge the latter may be at a lower effective 
temperature near the plate than near the cathode. 
Whatever the cause of this discrepancy may be 
the experiments prove with certainty that the 
thermal noise in these tubes is even lower than 
that calculated by Eq. (2). 


DISCUSSION OF RESULTS 


Moullin and Ellis? have recently contended 
that neither of the effects exists which have here 
been described and measured. They believe that 
the flow of electrons is random even in the 
presence of space charge, so that space charge 
should have no effect on shot noise, the reduction 
in noise at high space currents being caused 
solely by the decrease in internal plate im- 
pedance of the tube; and that there is no possibil- 
ity of the electrons in the space current sharing 
thermal energy with molecules, so that the in- 
ternal impedance of the tube cannot produce 
thermal noise. According to the views of Moullin 
and Ellis the factor j(07/d7)* in Eq. (1) should be 
replaced by 7, the other factors remaining the 
same. 

The dashed curve in Fig. 2 gives the noise 
calculated in accordance with this relation which 
denies the existence of the cushioning effect of 
the space charge. The agreement of the experi- 
mental points with the solid curve calculated 
from Eq. (1) shows that space charge is effective 
in reducing shot noise. The experiments prove, 
further, that the impedance of the tube con- 
tributes to the total thermal noise of the circuit 
but that in place of the cathode temperature an 
even lower temperature must be used. This tem- 
perature appears to lie about midway between 
the temperature of the cathode and room 
temperature. 

In conclusion I wish to express my thanks to 
Dr. J. B. Johnson for the helpful discussion and 
criticism given during the course of this work. 


7E. B. Moullin and H. D. Ellis, I.E.E.J. 74, 323 (1934). 
F. B. Llewellyn (1.E.E.J. 75, 395 (1934)) has discussed 
this article. 
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The Potential Distribution Between Parallel Plates and Concentric Cylinders Due to 
Any Arbitrary Distribution of Space Charge 


W. G. Capy, Wesleyan University 
(Received October 29, 1934) ’ 


Equations are derived for the electric potential and field intensity due to a distribution of 
space charge in a slab of dielectric between two infinite parallel conducting plates. The charge 
density may be any integrable function of the coordinate perpendicular to the plates. The 
corresponding solution for concentric cylinders is also given. A number of special cases are 
discussed, including that of a uniformly distributed charge, and one of sinusoidal distribution. 
The latter has an application in the theory of the piezoelectric resonator. 





N the author’s recent development of the 
theory of the piezoelectric resonator, which 
will soon be published, it has been found neces- 
sary to take into account the reaction of polariza- 
tion charges in the crystal upon the performance 
of the resonator. To this end, general equations 
have been derived for potential and field intensity 
due to a distribution of space charge in a slab of 
dielectric between two infinite parallel conducting 
plates, the charge density being any arbitrary 
integrable function of the coordinate perpendic- 
ular to the plates. This solution is then applied 
to the resonator problem. The derivation is given 
in the present paper, in the hope that the results 
may be useful in a wider field. Though it has no 
present bearing on the resonator problem, the 
corresponding theorem for space charge between 
concentric cylinders is also included. Special cases 
of the general problem in its relation to thermi- 
onic phenomena and discharge in gases have been 
discussed by many writers.' 


I. PARALLEL PLATES 


In Fig. 1 the x-axis is perpendicular to the 
plates, which are located at the origin O, and 
at X. The slab of dielectric, which has a dielectric 
constant k, occupies the space from A(x=a) to 


1 For example, C. D. Child, Phys. Rev. 32, 492 (1911); I. 
Langmuir, Phys. Rev. 2, 450 (1913), also 24, 49 (1924); W. 
Schottky, Verh. d. D. Phys. Ges. 16, 490 (1914), also 
— d. Radioakt. u. Elektronik 12, 147 (1915); T. C. 

ry, Phys. Rev. 17, 441 (1921); R. Gans, Ann. d. Physik 
69, 385 (1922); I. Langmuir and K. T. Compton, Rev. 
Mod. Phys. 3, 191 (1931); A. Trost, Phys. Zeits. 35, 725 
(1934); S. Werner, Zeits. f. Physik 90, 398 (1934); H. J. 
Van der Bijl, The Thermionic Vacuum Tube, New York 
1920, pp. 54-57; J. J. Thomson, Conduction of Electricity 
Through Gases, 3rd ed., Vol. 1, 1928, p. 370 ff; Wien- 
Harms, Handb. d. Experimentalphysik, 1928, Vol. 13, 
Part II, pp. 244 ff. 





B(x=a+e). On each side of the dielectric is a 
gap having k=1, devoid of space charge, the gap 
lengths being a and b, respectively. The space 
charge in the dielectric, represented by the curve 
p in Fig. 1, may be any arbitrary function of x, 
but is uniform in all directions perpendicular to 
the x-axis. The potential is assumed to be zero 
at O and X. If there is a potential difference 
between the plates, it can be superposed on the 
space charge solution. For a one-dimensional 
problem, Poisson’s equation has the form 


d’ V(x) /dx*= —(4x/k) p(x), (1) 


where V(x) and p(x) are functions of x alone. 

The following notation will be used, wherein 
C, and C are constants to be determined from 
boundary conditions. 


J o(x)ds= p(x) +Cs, (2) 
J ex(xdse= pals) +Cs, (3) 
s=k(a+b)+e, (4) 


where s may be called the “electrical distance” 
between the plates. Upon integrating Eq. (1) we 
have 


d V(x) 4r . 
——— = ——|pi(x) + Ci} = — E(x), 
dx k 
[a<x<(a+e)], (5) 
4n 
V(x) = — a beste) + Cit Cal, 
[a<x<(ate)]. (6) 
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Fic. 1. Space charge distribution in dielectric between 
parallel plates, with gaps a and b. 


These expressions for potential V and field in- 
tensity E are valid only in the region occupied by 
space charge, i.e., from A to B in Fig. 1. The 
values in the gaps a and 3 are easily derived as 
soon as the potentials at A and B have been 
determined. 

C, and C2 are found from the boundary condi- 
tions at A and B. At A, where x=<a, the electric 
intensity in the gap is — V(a)/a, while in the di- 
electric it is — V(a)/ka. In passing, it may be 
remarked that this recognition of the change in 
intensity at the surface of the dielectric obviates 
the necessity of introducing explicitly the 
polarization surface charge. Upon setting x=a 
in (5), we have 


dV(a) V(a) 
dx ka 





4r 
owe ted. 


From this equation, together with a similar one 
for the point x=(a+e), and the two equations 
obtained by substituting a and (a+e) re- 
spectively for x in Eq. (6), we can solve for Ci, 
C2, V(a), and V(a+e). The expressions for C, 
and C2 are 


1 
Ci=~—|{¢2(a) —kapi(a) — p2(a+e) 
Ss 
—kbpi(at+e)}, (7) 
1 
C2=-[(a+e+hkb) {kapi(a) — p2(a)} 
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—a(k—1){kbpis(a+e)+p2(a+e)}]. (8) 


If, in addition to the considerations discussed, 
the plate at X, Fig. 1, is at a potential Vi, we 
simply add the corresponding term to Eq. (6). 


POTENTIAL DISTRIBUTION AND SPACE CHARGE 


II. CONCENTRIC CYLINDERS 


In Fig. 2 the origin is at O, the axis of the 
cylinders. The two cylindrical electrodes, of in- 
finite length, both at zero potential, have radii 7; 
and 74. The space charge, distributed uniformly 
around and in the direction of the axis, is com- 
prised between the clyindrical surfaces of radii 
re and r3. It may be expressed as any arbitrary 
integrable function p(r) of the radius. The dielec- 
tric constant is assumed to be unity throughout 
the system. 

Poisson’s equation for the cylindrical case is 


id tk 





-—(r 


rdr dr none ©) 


Upon integrating and introducing the notation 


[roar =F\(r)+Ci, 


far F\(r)dr= F.(r) + Co, 


we have for the intensity 


dV(r) 4n 
— E(r)=-——= ——[Fi(r)+Ci} 
d r 


’ [re<r<r3], (10) 
and for the potential 
V(r) = —4r{ Fo(r) + Ci log r+ Co} 
[re<r<r3]. (11) 


The constants are determined by applying to the 
regions 772 and 73f4, which are devoid of charge, 
the usual equation for potential difference be- 
tween infinite concentric cylinders. Remembering 
that V(r) = V(r4) =0, we have V(re) = —reE(re) 





", 


Fic. 2. Partial cross section of concentric cylinders of 
radii 7; and 74, with space charge occupying the region from 
%o tO 73. 
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log r2/r;, and a similar equation for V(r3). From 
these equations, and by substituting 72 and r; for 
rin (10) and (11), we find 


C,\=— neal F.(re) gs F.(r3) 
log (r4/r:) | 


Yo Ys 
—log —F,(re) —log —Fi(rz)}, (12) 


r ’3 


C2= log (r2/ri) Fi(re) — Fo(re) — C, log ri. (13) 


EXAMPLES 
(1) Parallel plates, gaps a= b=0, R= 1 


The space charge may have any distribution 
p(x) perpendicular to the plates. From Eqs. (6), 
(7), and (8), the potential at any point is 


V(x) = (42/e) {xp2(e) +(e —x)p2(0) —epe(x)}, (14) 


where e is now the distance between the plates. 
In the special case where p is uniform this 
reduces to the usual parabolic equation 


V(x) = 2apx(e—x). 


(15) 


If the space charge has the form p(x) =po+cx 
the solution is 


V(x) = 3 ax(e—x) {3po+c(e+x)}. (16) 


This case is illustrated in Fig. 3, which is 
drawn on the assumption that po=0.5 e.s.u., 
e=1 cm, and c=1. The maximum value of V 
comes at x=0.54 cm. 
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Fic. 3. 


Linear distribution of space charge between 
parallel plates, with no gaps. 
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(2) Parallel plates, sinusoidal distribution of 
‘space charge 

This case is encountered in the piezoelectric 
resonator. For illustrative purposes it suffices to 
assume that we have a plate of piezoelectric 
crystal AB, Fig. 4, of thickness e and dielectric 
constant k, of relatively great extent in directions 
perpendicular to the X-axis, and with gaps of 
lengths a and b between it and the electrodes, 
which are at A’ and B’. The origin is now at O, 
Fig. 4, in the center of the dielectric. Owing to its 
state of instantaneous strain while in vibration 
in the X-direction, the plate is polarized parallel 
to X according to the equation P(x) = Py cos rx/e. 
This polarization, represented by the curve P 
in Fig. 4, gives rise to a space charge p of the form 
p(x) = —dP(x)/dx=po sin rx/e, where po =7Po/e. 
The quantities p:(x) and pe(x) in Eqs. (2) and (3) 
are now written 

wx poe? WX 


poe : 
pi(x)= —— cos—, p2(x)= ——sin —. 
T e 7 e 


(17) 


By the use of these expressions, and with due 
regard to the shift in the position of the origin, 
Eqs. (7) and (8) become C,=2poe?/m?s and C2 
= —kpoe?(b—a)/x’s. The potential equation is 
therefore, from (6), 


A poe” mx 2x k(b—a)) 


sin ———+ : 
1k e Ss Ss J 


[a<x<(a+e)]. (18) 


As is shown by this equation and represented 
in Fig. 4, the potential function has the form of a 
sine curve superposed upon a straight line. 

The field intensity in the dielectric has the 
equation 


E(x)= — 








V(x)= 


mx 2e 
cos ——— 
é TS 








d V(x) Apoe 
= ( (19) 
dx k 
The intensity is, under the conditions shown in 
Fig. 4, negative except near the bounding faces 
of the dielectric, where it becomes positive. It 
passes through zero at points corresponding to 
maximum and minimum potential. These points 
+x, are given by the equation cos 7x,,/e=2e/ms. 
In Fig. 4 lines are also drawn representing the 
values of E and V in the gaps a and b. If a=), 
the potential curve becomes symmetrical, passing 
through the origin. 





POTENTIAL 























Fic. 4. Sinusoidal distribution of space charge in dielectric 
between parallel plates. 


(3) Concentric cylinders, no region devoid of 


space charge 


In Fig. 2, a=b=0, 11 =7e, r3=75. The p distri- 
bution extends from 7 to rs. By making these 
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substitutions in (11), (12) and (13), we find for 
the potential at any point 


4nr 
V(r)= Fo(r4) — Fe(ri l 
(r) ee (rs) (r1)} log r 


rT 
—log —F2(r) —log r1F2(r.) +log reFalrd) | (20) 


el 


If, under the same conditions, p is uniform be- 
tween the cylinders, the functions F,; and F¢2 
become F,(r)=pr?/2, Fe(r)=pr?/4. Hence (20) 
becomes 


mp(ra—r;*) 
V(r) =—————— log r— apr’ 
log (74/11) 


mp(ri? log rs—r2 log 171). 





(21) 
log (74/11) 

This equation has the same form as Eq. (217) 

of Langmuir and Compton,! and in addition it 

gives the values of their coefficients A and C. 
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Potential Distributions About an Infinitely Extended Line Electrode on the Surface 
of a Horizontally Stratified Earth. II 


Morris Muskat, Gulf Research and Development Corporation, Pittsburgh 
(Received October 5, 1934) 


The analytical theory is given for the potentials and 
gradients to be found on the surface of a horizontally 
stratified earth in which the current source is an infinitely 
extended cable or line electrode lying on the surface. 
Power series expansions in the distance from the electrode 
are derived for small distances and asymptotic expansions 
for large distances. Detailed treatments are given for the 
two layer earth, for all thicknesses of the surface layer, and 
for the three layer earth in which the thickness of the second 
layer either equals or is twice as great as that of the surface 
layer, as well as formal expansions for general many layered 


earths. Simple closed hyperbolic function expressions are 
derived for the surface gradients for three layer systems in 
which the lowest layer is either of infinite or zero conduc- 
tivity. The solution for the direct interpretation problem in 
which the continuous vertical conductivity variation may 
be determined from a knowledge of the surface gradients is 
briefly outlined, following the method of Slichter and Langer 
for the case where the current source is a point electrode. 
The analytical results of the indirect interpretation problem 
for a stratified earth are illustrated by thirty numerical 
and graphical examples. 





INTRODUCTION 


HE problem of electrical prospecting by 

means of point electrodes has recently been 
the subject of considerable study both from the 
point of view of indirect! and direct? interpreta- 
tion. While much remains yet to be done in the 
development of convenient and rapid numerical 
techniques, especially in the case of the direct 
interpretation method, the analytical theory of 
the problem may be considered to be fairly well 
completed insofar as one assumes a uni-dimen- 
sional variation of the conductivity. However, 
the case of the analogous two dimensional 
system in which the current or potential source 
is an infinitely extended cable, and which has 
been frequently used in ore prospecting, has not 
been given any detailed treatment. For the 
stratified earth, several cases have been discussed 
graphically by means of series of logarithmic 
terms* while integral expressions have been 
derived by King! for the surface gradients for 
continuous variations of conductivity which are 
either exponential or parabolic. 

Although, as King points out, the theory of 
the line electrode can be formally derived from 
that for a point electrode by integration of the 
results of the latter and will give no essentially 

'L. V. King, Proc. Roy. Soc. A139, 237 (1933), and Phil. 


Trans. Roy. Soc. A233, 327 (1934); M. Muskat, Physics 4, 
129 (1933). This latter paper will be referred to in the text 


as I, 

? L. B. Slichter, Physics 4, 307 (1933); A. F. Stevenson, 
Physics 5, 114 (1934), 

*L. J. Peters and J. Bardeen, U. of Wisconsin Engineer- 
ing Expt. Station Bul. 71, 1930. 
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new information beyond that of the point 
electrode,‘ the line electrode gradients may be 
expressed as series of exponentials rather than 
Hankel functions and hence are better suited for 
practical interpretation when one of the layers 
is of zero or infinite conductivity. In fact, King 
indicates how to derive from empirical point 
electrode gradients the corresponding line elec- 
trode gradients which are more convenient for 
analysis in terms of the electrical depth con- 
stants. 

While most of the analysis will be ultimately 
expressed here also in terms of the surface 
gradients, the potentials themselves will diverge 
if obtained by simple integration of the integral 
expressions for the point electrode potentials. 
Furthermore, when none of the layers is of zero 
or infinite conductivity the point electrode 
potentials or gradients cannot be expressed as 
Hankel function series and power series must be 
used for small distances from the electrode; 
these will evidently diverge on integration. We 
shall, therefore, briefly outline the direct deriva- 
tion of the integral expressions for the potentials 
of the line electrode system in a two layer earth 
without the use of the point electrode expressions. 
The expressions for the 3 layer earth will be 

4On the other hand it should be observed that the two 
dimensional prospecting problem implies no physical 
simplification as compared to the point electrode system, 
as in both cases the earth is supposed to be horizontally 
stratified with uniform conductivities in planes parallel to 
the earth’s surface. The reduction of the analytical problem 
to one of two dimensions is based on the use of a theo- 


retically infinite cable as the current source rather than a 
point electrode. 
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given without derivation and only for the 
surface gradients, and may be considered as 
obtained by integration of the point electrode 
gradients or by independent solutions similar to 
that given for the two layer earth. Finally, it 
may be mentioned that although the case of 
continuous conductivity variation appears to be 


more general than that of‘ discrete layers of 
different conductivity, the latter representation 
probably corresponds more closely to the situa- 
tions occurring in practice and cannot be 
obtained from the continuous conductivity vari- 
ation theory, but must be derived by a separate 
analysis. This analysis will be given here. 


I. A Two LAYERED EARTH 


The analytical problem may be formulated 
as follows: find potential functions #; and 
pertaining to regions (1) and (2) (cf. Fig. 1), 
such that: 


= log {x*+(y—h)*} :lim {x?+(y—h)?} 0; 
db;/dy=0 : y=0; (1) 


Di= $25 019%1/dy=o29s/dy : y=}, 


where the unit of length has been taken as twice 
the thickness of the upper layer and «a, o2 denote 
the conductivities of the two layers. 





. b(0, h) 
/e (1) G, 


| 





(2) 





| 


y 


Fic. 1. Diagrammatic representation of a line electrode at 
(0, 2‘ in a two layer earth. 


Although it is not difficult to sum the series of logarithmic terms corresponding to the system of 


images giving the solution to the Eqs. (1) the direct integral method has the advantage that it 
may be readily generalized to systems containing any number of parallel layers, and even to those 
in which the conductivity varies continuously. Here one may begin with an integral representation 
of the primary source which may be written as: 


co da e72(h-ya rV< h 
$= log {x*+(y—h)*} =2f | ete—cos dna || 
0 @ ghee s 935 


(2) 


The perturbation in the potential in the upper medium due to its closure at y=0 and the presence 
of the lower medium at y=} may then be expressed as: 


© da 
,= 2 f —[Wila)+ {A e244 + Bierve} cos 2xa | (3) 
0 a 


and the total potential in the lower medium as: 


P.=2 f (da/a)[WV2(a)+A 2e~@u-D@ cos 2xa ], (4) 
9 


where VY; and W2 are introduced to take care of the term e-2* in @9. Applying the boundary conditions 
of Eqs. (1) to $;,=&)+%, and and adjusting VW; and W2 so as to insure the convergence of Eqs. 
(3) and (4) it is found that: 


AA 1= 2(0;/02) cosh (1 —2h)a+2 sinh (1 —2h)a; 
AA o= 401/02 cosh 2ha; 
Wi(a) = (20;/02—1)e-**; 


AB,= —2(1 —a;/o02)e-* cosh 2ha, 
A= —2{cosh a+o;/e2 sinh a}, 


Wo2(a) = 201/a2-e7°*, 


(5) 





——— 
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so that denoting by k the quantity: (0,—2)/(oi1+ 02), 


o;/02<1 :k=—e™?: 


© daft cos 2xa 
a=4 f san pen ye 24 — —_—_—_—_ 
0 


sinh (a+7—2ha) cosh 2ya : | 























a cosh (a+) sinh (a+y—2ya) cosh 2ha : y>h 6) 
© daft sinh y cosh 2hae~@4-¢ cos 2xa r 
?,=4 f —| tanh ye?*+— | 
0 @ cosh (a+v7) 
o;/02>1:k=e 7: 
© dat cos 2xa cosh (a+7—2ha) cosh 2ya : y<h) 7) 
a=4 | —{ ctnh ye?*— | 
Jo al sinh (a+~) cosh (a+7—2ya) cosh 2ha : y> 
L 7) 
© daft cosh y cosh 2hae~@4-) cos 2xa 
nas f —j ctnh ye **— . =| 
._ 2 j sinh (a+7) 





Although Eqs, (6) and (7) remain valid when h=0° it is convenient for the practical discussion 
of the potential distributions to first separate out ®p. Restricting, for the sake of brevity, the further 
treatment to this case of most practical interest where 4=0 and considering, furthermore, only the 
surface potentials #;(y=0), Eqs. (6) and (7) may be transformed to: 























e~? po e- “day cosh y cos 2xa 
o1/o2<1:k=—e 7: ,(y=0)=log x— f : er — | (8) 
cosh y 79 a L cosh (a+v7) 
e~? e- “day sinh y cos 2xa 
o1/o2>1:k=e7: (y= 0) =log x+— f ge——— | (9) 
sinh y “9 a L sinh (a+y7) 
where a factor of 4 has been taken out. oi/oz<1; k= —e*7; 
Due to the divergence of the separate parts ” @-* oie Son de 
of the integrals in Eqs. (8) and (9) they are &’=1/x—2e-7 f ; (10) 
unsuited for detailed discussion. As indicated in 0 cosh (a+) 
the introduction we shall, therefore, remove the , — ee 
ri f 0;/d2> 1 P4 k=e°7; 
constant terms of the integrands by differentia- 
tion with respect to x; convergent integrals are ae _ ee sin 2xa da 
then obtained which are, in fact, equivalent to &'=1/x+2e% f ; - (11) 
0 sinh (a+y¥) 


Fourier integrals. Although the original problem 
as stated in Eq. (1) referred to the potential 
distribution itself, there is no real loss in gener- 
ality in expressing the solutions in terms of the 
gradients of the potential as the series to be 
derived for the gradients may be readily inte- 
grated term by term. Differentiating, then, Eqs. 
(8) and (9) and dropping the subscript in 4%, 
we have for our final integral expressions for 
the surface gradients: 


5 The apparent non-vanishing of d%,/dy at y=0 for h=0 
is due to the discontinuity of the integral representation of 
0% 0/dy (cf. Eq. (2)) and hence of d4,/dy at y=0. When 
is first separated out the difficulty disappears. 


A. Expansions for small values of x 


By a method analogous to that developed in J, 
the expansions of Eqs. (10) and (11) for small 
values of x may be readily derived by expanding 
first sinh (a+) and cosh (a+) in series of 
exponentials, then substituting the power series 
for sin 2xa, and finally integrating term by 
term. Thus it is found that both for o;/o2 > and 
€i: 


= 1/x+2¥ (—1)x2™+1 DS k/n™42, (12) 
0 1 
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Here again as in J, the radius of convergence 
of the series may be advanced from the value 1, 
by summing out the terms with n=1, 2, etc. 
Doing this for the first two, &’ may be rewritten 
as: 


1 2kx 2k?x 


J 425° (—1) "fn (b)x2"41, (13) 
4+x° 0 





+ + 
x 14+ 
where f»,(k)= Sk" /nemt?, 

3 
B. Expansions for large values of x 


By direct application of the general formula 
(cf. I Eq. (46)): 
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Need at (—1)™f?"(0) 
f f(a) sin xa da= >> —————-__ (14) 
0 0 xemtl 


it follows at once that the asymptotic expansions 
of the surface gradients are: 





01 loa< 1 . 
- (- 1)” gem 
@’=2>)° (tanh y); (15) 
0 (22)°"** oy" 
o1/02> es 
a (— 1)” ozm 
p’ = (ctnh y). (16) 





0 (2x)?m+1 oy?" 


Carrying out the differentiations for the first 
four terms it is found that for both o;/e2 > and 
<i, 


01 1 1 if 1 
=| 14 ast/ort) — (0 —ast/o2) (Bot of —2)+ (1 o4/02)(17 ~ 000 8/a2+ 4501/00 
2x? 2x4 Ax6 


O2X 


1 


It is of interest to note that except for the 
numerical factors the variation of ®’ with o/c 
is exactly the same as that for the corresponding 
point electrode system (cf. J, Eq. (16)). And as 
in that problem the asymptotic expansions fail 
both for o.= © and o2=0, due to the degenera- 
tion of the inverse power variation of ®’ into an 
exponential variation. 

Thus for o.=*, y=0 for Eq. (10) and ®’ 
becomes : 


© e-* sin 2xa 
w= 1/x-2 [ ——da=rcsch rx, 
0 cosh a 


(18) 


which evidently behaves as e~** for large x. 
Similarly for o.=0, y=0 in Eq. (11) and 


© e-*sin 2xa da 
~’=1 x2 f =nctnh rx, (19) 
0 sinh a 





which approaches 7(1+2e-°**) for large x. The 
constant term here, which represents a linear 
term for @, corresponds to the logarithmic term 
in the asymptotic expansion for the point 
electrode problem (cf. J, Eq. (18)). 

The gradient distributions for o:/02=0, 0.20, 


—e —o1"/027)(31501°/ 02° — 3600 1'/o2'+3780;?/o22—62)+--- | (17) 
x 





1.0, 5.0 and ~ are plotted in Fig. 2, where for 
convenience x®’ has been used as the ordinate 
rather than ®’ itself. For small x Eq. (13) was 
found quite satisfactory as was Eq. (15) for 
large x for o:/02=+4. However, for o:/o02=5, Eq. 
(16) was first modified to the form: 


o;/o2>1: 
©’=2 sin 2xyCi(2xy) —2 cos 2xy{ Si(2xy) — 2/2} 
co (—1)™ | am 


m 
0 (2x)?™+1 


T'(2m+1) 
(ctnh y) — 


2 sithdinendanidiniiiiaaiiii 
oy" ym 





|, (20) 


which is analogous to J (15a). The convergence 
was then so rapid that with only 3 terms of the 
series the remainders were less than 0.003 for 
x25. 

The curves here for x’ are strikingly similar 
to those for x® in the case of the point electrode 


¢ (cf. I, Fig. 3). The only differences are that here 


the approach to the asymptotic values o;/02 is 
somewhat slower than in J and that the slopes 
of the x®’ curves vanish at the origin whereas 
those of x® for the point electrode system have 
non-vanishing slopes at x=0 and hence show no 
inflections. 
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x 


Fic. 2. Surface gradients (#’) in a two layer earth. 
o,=conductivity of surface layer; o:=conductivity of 
lower layer. x = (horizontal distance from line electrode) / 
(twice the thickness of the surface layer). 


II. AN EARTH WITH MorE THAN Two LAYERS 


A. General theory 

For a general many layered earth, the direct 
integral method leads to an expression for the 
surface gradient of the form: 


£'(x)= 1/x+2f f(a) sin 2xada, (21) 
0 


where f(a)~e~** for large a, thus insuring the 
convergence of the integral. For small values of 
x a power series in x may be derived for &’—1/x 
by expanding sin 2xa in a power series. Thus it 
is found that: 


o> (—1)™(2x)?m+1 pm 


f a"+1f(a)da; (22) 
T(2m+2) Jo 





1 
©’ (x)=-+2 
x 0 
the integrals still remaining may be evaluated 
by first expanding f(a) in a series of exponentials 
and then integrating term by term. 
For large values of x, a partial fraction ex- 
pansion of f(a) may be used. Thus f(a) may be 
expressed as: 


ectrode : 





(1) 0, 














Vy 


Fic. 3. Diagrammatic representation of a three layer earth. 


f(a)=f(0)+Lbn{1/(a—an)+1/an}, (23)° 


where a, are the poles (simple) of f(a) and dD, 
the corresponding residues. ®’ then takes the 
form: 


®'(x)=1/x{1+f(0)+20b,/an} 
+>0,G,( 7 2xan), 


sin xa da 
eg 
0 atan 


To get an asymptotic expansion for ®’ it 
follows directly from Eq. (21), on applying 
Eq. (14), or from Eq. (24) on observing that the 
asymptotic expansion for G is: 

—— 1)"1(2m+1) 
ry (z)2m+1 


wo (—1)"f?"(0 
at /e47F nf (0) 


9 (2%)**t 


Furthermore, when the a, are pure imaginaries 
the functions G satisfy the relation: 


G(ixan) +G(—ixa,)=27re~*, (28) 
so that if b(ta,) =b(—ia,), Eq. (24) becomes: 
(x) =1/x(1+f(0))+27d0b,e-%7%". (29) 


6 When a=0 isa pole of f(a) the term bo/a must first be 
subtracted from f(a) before the expansion. 

7Use has been made here of the apparently incorrect 
result that fo sin 2xa da=1/2x. The difficulty is only 
apparent, however, and arises from the integral representa- 
tion of Eq. (2) which leads to the same integral and value 
for y=h, on differentiation of #9. A convergence factor as 
e~** in Eq. (21) before the expansion and integration would 
evidently lead to Eq. (24) after making ¢« approach zero. 


(24)? 


where (25) 





(26) 





that (27) 
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B. The three layer earth 


An analysis similar to that given in Part I, shows that for a three layer earth (cf. Fig. (3)) f(a) 
has the form: 


2[ koe~**-+ kse 2 ] 
f(a)= (30) 
1 — koe~** + k3e?2"( —1 + kee**) 





where 
ko= (o1—02)/(o1+02); k3= (¢2—03)/(o2+03); h=he/hy, (31) 


o1, 92, 03 being the conductivities of the top, intermediate and lowest (infinite) layers, 4; the depth 
to the second layer, hz the depth to the third layer and a@ has been so chosen that the unit of length 
for x in Eq. (21) is 2 dy. As the f(a) here is exactly the same as that for the point electrode problem 
(cf. I (25)), the results of the latter, which depend only on the form of f(a), may be immediately 
applied here. 


Thus for the. first two terms in the asymptotic expansion of Eq. (27) for general values of 1, 02, o3 
and h we get: 


xb’ (x) = = i-— 
O3u 202703°x? 





a} (o2"—a3") {o1—o2+o2h}?+037(0; — 02") 
re ‘| (32) 
which is identical with the corresponding expression of J(26) except for a factor of 2 in the second 
term. It follows from Eq. (32) that for very large x, x’ approaches the value o;/o3 and is independent 
of oe. If o1=03, xb’ will, therefore, assume the same value for large x as in a homogeneous earth. 
However, its approach to the limiting value 1 will depend on oe, and will be a descent from higher 
values if o2<o; and rise from lower values if o2> 01. 
For small values of x the power series expansion of Eq. (22) may be used to compute ®’(x). How- 


ever, in the evaluation of the integrals /°a?"*'f(a)da it is convenient to treat separately the cases 
corresponding to different values of /, as will be seen presently. 


Power series expansions for special values of h 


h=2. When the thickness of the second layer is exactly equal to that of the surface layer (h=2), 
f(a) can be expressed as: 


fla)= 





k3e-3*+ koe~* e711 e-v2 
| (33) 


[ks*(1—ks)?+4ks]}!|sinh (a+7;) sinh (@+72) 
where 2k3e27%1= — ko(1—k3)+[ko?(1 —ks)?+4ks |? =2k3/s1, (34) 
2k3e272= —k,(1 —k3) act [Rke?(1 —k3)?+4ks3 ]*=2k3/Sso. 


Expanding now Eq. (33) in exponentials and integrating Eq. (22) term by term it is found that: 








: 2 © roy ke ks 
B(x) =1/2+ F(—1) mS (54"—s9") ( + ). (35) 


[k?(1—ks)?+4k; |} 0 1 41)2m42 


the radius of convergence of which can be extended by summing out the terms with »=1, 2, etc. 
h=3. When the second layer is twice as thick as the surface layer, f(a) can be expressed as: 


s3e~71/(Se—S1)(S3—S1) $9¢~72/(S1—S2)(S3—S2) S$ 3€~7%8/(S1—S3)(S2—Ss) 


sinh (a+71) sinh (a+72) sinh (a+7s) 





fla)=e(br+ bxe-)| , (36) 


where s;=e~?”' and the denominator of f(a) has been factored in the form II(e~?*—e?”') so that the 
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7; or s; can be obtained by the solution of a cubic equation. Expanding f(a) as a series of exponentials 
and integrating Eq. (22) term by term, ®’(x) takes the form: 








2 co co ko ks 
©’ (x)=1/x+ > (—1) "2h! + 
(S2—5$1)(S3—51)(S3— Se) 0 1 n2mt2 (n+ 2)2+2 


+ {51"*1(S3— Se) +52"*1(51— 53) +53"*"(S2—51) }. (37) 


Here, too, a preliminary summation with respect to m of the terms for n=1, 2, --+ will extend the 
radius of convergence of the power series. 

In a similar manner power series expansions can be derived for systems in which h is any integer 
or rational fraction, although the degree of the equation to be solved for the s; would equal h, if 
it is integral or to the numerator of h if it is a rational fraction. Thus for h=3/2, there would again 
be three s; just as for h=3, which would have to be obtained from the roots of a cubic equation. 


Lowest layer of zero conductivity® 
When the lowest layer is of zero conductivity (¢;=0; ks=+1) it may be easily shown that all 
the poles of f(a) are imaginary and are in fact given by those values a, which satisfy the equations: 


[sinh a,h]/[sinh a,(h—2)]=ke (38) or sinh 2a,=0, and sinha,hk=0, (38a) 


the latter of which requires that h=2p/q, p and q being integers. . 
As a=0 is a pole of f(a) for ks= +1, the expansion of Eq. (23) must be modified to the form: 


1+k: 
Qa, k3=1)= —1 b, 1 —~An Ans, 
f(a, ks=1) (tht bbe + Lbal{l/(a— an) +1/an} (39) 


sinh 2en(h—1) 

















where b,=— for Eq. (38) 
sinh 2a,(h—1)—(hA—1) sinh 2a, 
40 
1+ kee? ( ) 
= — for Eq. (38a). 
h+(2—h)kse?= 
Hence, as b(ta,) =b( —ian) 
Oe in eS 
ee > ” izan. 41 
2koth—koh (41) 
where only the positive imaginary values of a, are to be included in the sum. 
Thus for h=2 the above equations give: 
a,=nri/2 :b,=(1+(—1)"ke)/2. 
Tv 
Hence &'(h=2; o3= o~~yum mx/2+ke tanh rx/2}. (42) 
For h=3: 
1+k» 1 
a,=nt1; b= : a,= +— cos (ke—1)/2+n77; b,=(1—k:2)/(3—ke) 
3—ke 2 
Tv 
and ©’ (h=3, o3=0)= {(1+:) ctnh rx+2(1—k:2) cosh (x—86)x csch rx}, (43) 


—Re 
where cos 0= (k2—1)/2. 


8 It should be noted that while the results presented zero or infinite conductivity is of finite thickness. For 
here for lowest layers of zero and infinite conductivity are such strata shield those above them from those below, 
derived from expressions referring to a three layer earth, and regardless of the stratification of the medium below 
thus implying that these lowest layers are of infinite these strata, the surface potentials and gradients will be 
thickness, they remain valid even when the third layer of | the same as if they were of infinite thickness. 
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For h=4, Eqs. (38), (38a) and (40) give: 





1+(—1)"ke 1 1—k,? 
a,=nni/2; 2= 3 an= +-— cos ko/2+n71; b,.= ’ 
2(2—(—1)"ke) 2 4—k,? 
so that 
m(1+ke) a(1 — ke) 27r(1 —k,’) 
©’ (h=4, o;=0)= h csch rx+ cosh (r—@)x csch rx, (44) 








—— ctnh rx +———_ 
2(2—ke) 2(2+k:) 4—k? 
where @=cos™ ke/2. 


It may be observed that the expressions of Eqs. (42), (43) and (44) reduce to that of Eq. (19) for 
the corresponding two layer earth when kz is set equal to 0 or 1. When k:= —1, they reduce to 
Eq. (18) for a 2 layer earth with a lower layer of infinite conductivity. 


Lowest layer of infinite conductivity 


When the lowest layer is of infinite conductivity (¢;= ©, ks=—1) the poles of f(a) are again 
imaginary, being given by the equations: 

[cosh a,h]/[cosh a,(h—2)]=ke (45) 

or sinh 2a,=0, and cosh a,h=0. (45a) 


Here, a=0 is no longer a pole of f(a@)—except in the trivial case ke=1, which reduces the problem 
to a two layer earth with the lower layer of zero conductivity—, and f(a) can be expanded as: 








f(a)= -14+ 200, {1/(a—an)+1/an}, (46) 
sinh 2a,(h—1) 1+Roe-2an 
where 0,= , for Eq. (45) = for Eq. (45a). (47) 
sinh 2a,(h—1)+(h—1) sinh 2a, h+(2—h)k2e?= 
As here, too, b(ia,)=b(—ia,), we get: &’(x)= 27> b,e7*72", (48) 


where the sum extends only over the positive imaginary values of a,: 
Applying these results to h=2 we find: 














1 
n= >. cos"! ke+nni; b,=}, 
so that ©’ (h=2, o3= ~)=7 cosh (r—8)x csch rx; cos 0= ke. (49) 
1 1+ke 1—ke 
For h=3: a,=+— cos (1+ke)/2+n71; b,= . an=(2n+1)7i/2; b,= 
2 3+ke 3+k: 
a(1—ke) 2r(1+hke) 
and &'(h=3, o3= © ) =————— csch rx -+————_ cosh (2 —6)x csch rx, (50) 
3+ke 3+ke 
where cos 6= (1+:2)/2. ; 
1 ko+(ko?+8)! 3ko+(k?+8)! 
Finally, for h=4: a,=+- cos{ |-+nxi :b,= 
2 4 4(k?+8)} 
—3ke+(k?+8)! 








1 ra 


|-+nxi ->,.= 


4 


4(k.?+8)! 
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and 
ma csch rx 

’’(h=4, ¢3™ 0 ) = ——_—_—_—_} (3ko+ (Ro? +8)!) cosh (3 —0:)x-+(—3ko+(keo?+8)?) cosh (31 —42)x}, (51) 
2(k2?+8)! 








“| 


| 
P ; 


4 


where 6,;=cos™! 


. 6.= cos} 


Here, too, the Eqs. (49), (50) and (51) reduce to the appropriate 2 layer earth expressions’ of 
Eqs. (18) and (19). For ke=1, they reduce to the value of Eq. (19) for a 2 layer earth with a lower 
layer of zero conductivity and for k:=0, —1 to that of Eq. (18), or its equivalent, for a 2 layer earth 
with a lower layer of infinitely high conductivity. 

In a similar manner the surface gradients for other values of # may be obtained by first solving 
Eqs. (38) and (45) for the a,, determining the b, by Eqs. (40) and (47) and finally ®’(x) by Eqs. 
(41) and (48). For rational values of / the series of Eqs. (41) and (48) will always be summable as 
groups of geometric series; when / is not rational the series will still converge rapidly even though 
they may not be summable in closed form. 

The same methods as used here can evidently be also extended to the case of systems with more 
than 3 layers. While it may be more difficult to prove, it will undoubtedly be found that the cor- 
responding f(a) will again have only imaginary singularities when one of the layers has a vanishing 
or infinite conductivity, and hence give rise to surface gradients which vanish exponentially if there 
is present a layer of infinite conductivity, or which exponentially approach constancy if one of the 
layers has a vanishing conductivity. Eq. (24) will be valid in any case, except for the correction 
indicated in the footnote when f(0)= «. 

The analytical results derived above are illustrated in Figs. 4, 5, 6, 7 and 8. In Fig. 4 the surface 
gradients multiplied by the distance from the current source are plotted against the latter for three 
layer earths for which o,=03 and o;/o2=}, 5, with /e/h,;=2. For small x, Eq. (35) was used with the 
terms for n=1, 2, 3, 4 summed out separately, s;, sz being real for og.=50, and complex for o,=5oze. 
For large x the extension of Eq. (32) was quite satisfactory for o;=5e2. In the case of ¢2=501, how- 


x 


Fic. 4. Surface gradients for the case of a slab of thickness equal to that of the top layer, immersed in an 
otherwise homogeneous earth. Notation same as in Fig. 3. 
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Fic. 5. Surface gradients for a three layer earth in which the lowest layer has a zero conductivity and the 
conductivity ratios of the two upper layers are 1: 2 and 2: 1. h,=thickness of surface layer; ho—h; 
= thickness of second layer; x = (horizontal distance from electrode) /2h,. 


ever, the analog of Eq. (20): 
ko(1—1/s;) 
~ [hP(1—hs)?+-4ks} 


U 





is 
[> 


m 
0 (23)" oy?" 


(—1)™ oem 1 
{etnh n-—| +2 sin 2xyiC1(2xy1) 
Y1 


ko(1 —1/se) wo (—1)"™ 9 
+008 2y3(n/2—Si(2ey,))} |- 3 (ctnh y2) (52) 

[ko?(1 —ks)?+4k; |! 0 (2x)? dy? 
was found very satisfactory, 51, se, y1 and yz being defined by Eqs. (34). Curves of x’ for 3 layer 
earths with lowest layer of zero or infinite conductivity and h=2, 3, 4 and o:/o2.=+}, 3, 2, 5, are 
given in Figs. 5-8. Those for o;= © are very similar to the corresponding ones for x@ for the point 
electrode (cf. Figs. 8 and 9 in J) except that here the x®’ curves have horizontal tangents at x=0. 

In the case of ¢;=0 the values of x®’ are all positive and approach the asymptotic forms: 





TX 
x’(¢3;=0)> 4 h=he/hy, (53) 
1+(h—1)o2/o1 





while the ® for the point electrode decrease monotonically and approach —~ logarithmically. 


Here, too, as will also be true in the general case,—as follows from Eq. (21)—the curves for x®’ 
have vanishing slopes for x=0. 


III. Direct INTERPRETATION PROBLEM FOR CONTINUOUSLY VARYING CONDUCTIVITIES 


To complete the discussion of the problem of the infinitely extended line electrode we shall briefly 
outline, following Slichter’s? analysis for the point electrode, the solution to the direct interpretation 
problem. Thus for the general continuously varying conductivity o, and uniform horizontal con- 
ductivity o,, the fundamental solutions for the potential distribution are of the form: 


cos 
P(x, y)=Vily,a)) . ax, (54) 


sin 
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Fic. 6. Surface gradients for a three layer earth in which the lowest layer has a zero conductivity and the 
conductivity ratios of the two upper layers are 1:5 and 5: 1. 


xo (x) 





% +o 
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Fic. 7. Surface gradients for a three layer earth in which the lowest layer has an infinite conductivity and 
the conductivity ratios of the two upper layers are 1 : 2 and 2: 1. 
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Fic. 8. Surface gradients for a three layer earth in which the lowest layer has an infinite conductivity and 
the conductivity ratios of the two upper layers are 1:5 and 5: 1. 


where Y; is the solution vanishing at infinity which satisfies the equation: 
Y,"+(¢,'/cy) Y\'—(o2/cy)e Y,=0. (55) 


Replacing now the cable by an infinite flat strip of width 2a at the center of which the potential 
is ®y, and which is the source of a current of strength C per unit length, the potential distribution 


will be: 
C in aof Y;(0, a) cosxaYV;,(y, 
OI 6 
rac(0)J) a L¥\0,a)  ¥1(0, a) 








The surface gradients are, therefore: 


sin aa Y,(0, a) sin xa 


#()=-— f : . = 
‘ ~ ae,(0) r Yi(0, a) ’ " 


which, being a Fourier integral, may be inverted as: 








aY,(0,a) —2aac, 


k(a)=-— —= few sin ax dx. (58 
Y,(0,a) Csinaady 





Thus here, the “‘kernel’’ k(«)—in Slichter’s notation—can be found by integrating the surface 
gradients ’(x) as a Fourier integral, whereas for the point electrode the determination of k(a) 
required a Bessel function integration over the surface potential itself. Knowing k(a), the con- 
ductivity o,(y) is then given by the Slichter-Langer theory in the form: 


oy(y) = aoe Sora ay, (59) 








where a;(y) is to be considered as formed by a Taylor expansion with coefficients a: (0) given by: 
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a,’"(0) 
2! 


a;(0)= a4; a,;'(0)= —a;?+2a2; a;>— 2a\a2+ 243, etc. 


and the a’s are to be obtained by expanding k(a) in the form: 
k(a)=1+4;/a+d2/a?+---. (60a) 


The writer is indebted to Mr. M. W. Meres for the numerical and graphical work of this paper 
and to the executives of the Gulf Research and Development Corporation for permission to publish it. 
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A method is given for solving logarithmic potential 
problems in which the potential is preassigned over part of 
the boundary and the normal derivative over the re- 
mainder. A “‘mixed” Green's function is defined such that it 
vanishes over part of the boundary of interest and its normal 
derivative vanishes over the remainder, and the solution 
for the logarithmic potential is expressed in terms of this 
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function and the boundary conditions by means of Green’s 
theorem. Taking the case of the infinite quadrant as the 
“‘prototype” region the solutions for other regions such as 
the infinite half plane, infinite and semi-infinite strips, and 
the circle are obtained by mapping them on the infinite 
quadrant; these are illustrated by several specific examples. 





LTHOUGH methods of treatment are 

available for the solution of mixed bound- 
ary value logarithmic potential problems in 
which the normal derivative vanishes over part 
of the boundary and the potential has constant 
values over the remaining parts of the boundary! 
no detailed discussion seems to have been given? 
of the corresponding problems in which the 
normal derivatives and potentials are pre- 
assigned arbitrarily. Thus even the simple case 
in which the potential is specified over the right 
half of the x axis and the normal derivative over 
the left half, does not seem to be solvable by the 
methods given in the standard texts. In the 
following, a method will be given for treating 


such problems, when restricted to two dimen- 
sions. 


The principle of the method is simply based 
upon a combination of the theories of Green’s 
functions and conjugate function transforma- 
tions both of which have found wide applications 


1 J. Jeans, Electricity and Magnetism, p. 352, 5th ed. 

? The writer is indebted to Professor H. Bateman for 
calling to his attention a series of papers on the present 
problem by Carl Neumann in the Lezpziger Berichte, Math. 
Phys. Klasse 61, 156 (1909). In the first paper of this series 
Neumann outlines two methods for the treatment of the 
mixed boundary value problem; one of these consists in the 
reduction of the general case to that in which dv/dn is a 
constant on I; the other is essentially identical with that 
given here, and Neumann sets up an equation equivalent 
to Eq. (3). However, he gives no applications or discussion 
of his equivalent of Eq. (3), although in the later papers of 
the series he develops his first method at great length after 
interpreting the reduced problem as that of finding the 
charge distribution on the segment IT, taken as a con- 
ductor, induced by an appropriately chosen external charge 
distribution. While several cases of the reduced problem 
are solved by Neumann by this latter method, the analysis 
required is not only much more complex than that invelved 
in the method developed here, but moreover, the process of 
“reducing”’ the original problem requires the preliminary 
solution of the ‘“‘Neumann”’ problem for the contour T 


with dv/dn having the value 0 on I; and its original value 
on Pe. 
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in the discussion of the Dirichlet and Neumann 
potential problems. The Green’s function used 
here will, however, be a ‘‘mixed”’ function rather 
than those employed heretofore. It is by the 
introduction of this type of function that the 
simple results obtainable by the method of 
conjugate functions alone! can be generalized, 
as will be seen presently. 

Thus from the Green’s theorems for two 
dimensions it is readily shown that the potential 
function v(x, y) within a contour [ may be 
expressed in terms of the boundary values of v 
and its normal derivative by the equation: 


dG(x, y; x’, v’) 





2rv(x, y)= flow y’) 
Jp 


on 
dv(x’, y’) 


on 


— G(x, y3 x’, y’) fs (1) 


where (x’, y’) represents a point on I, m is the 
exterior normal and G is a function of the point 
pair (x’, y’) and (x, y) and harmonic in each 
except at the point (x, y) =(x’, y’) where it has a 
logarithmic singularity. 

If now the countour I be separated into two 
parts I’, and [2 and G is chosen so that: 


G=0, on T;; 0G/dn=0, on To, (2) 
G may be defined as the “mixed’’ Green’s 
function for the contour I divided into T, and 


T;. Eq. (1) becomes then: 


1 0G 1 ov 
v(x, y=— f v—ds—— }| G—ds, (3)? 
2r~7ry On 2rv¥re On 


so that if G has been found, the value of v at 
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(x, y) may be computed from a knowledge of v 
on T, and of dv/dn on Tz. 

lr, and IT, may, of course, be themselves 
composed of several segments and Eq. (3) will 
still remain valid. Furthermore, the contour may 
be finite, or infinite and bound infinite regions. 
The rigorous proof of these features as well as 
further properties of G, such as its symmetry in 
(x, vy) and (x’, y’), may be carried through by 
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the methods already developed for the study of 
the ordinary Green’s functions, and will not be 
reproduced here. 

The problem that remains is the determination 
of G. This will be accomplished by first setting 
up G for simple regions and then treating more 
complicated systems by mapping them on the 
“prototype” regions by means of conjugate 
function transformations. 
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A. T; AND T, ARE SINGLE SEGMENTS 
(1) The infinite quadrant 


The prototype for regions in which I’; and I: are each composed of single segments will be taken 
as the infinite quadrant for which G can be found by a simple superposition of two source and two 
sink potentials. . 

Denoting the plane for the infinite quadrant as the ‘“‘w plane’”’ where w= +77 and taking for I), 
I, the positive — and 7 axes (cf. Fig. 1) we have here at once: 


t(€o— €)?+(mo— 0)?} {(Eo+ £)? + (n0— 0)*} 











G= > log , (4) 
((€o— €)?+ (mot n)*} | (Eo+£)?+(n0+7)*} 
so that: 
n 1 1 a cs 2+ (no— 1)? 
oe == f 160} + {at-— f° Fim) log dm, (5) 
7/0 (fo—£)*+? (Eo+£)? +7? 27 Vo 2+ (not+n)? 
where here and below f(£) and F(n) denote the values of v and dv/dn on T, and I>. 
If: F(o)=0; f( 0) =o, (6) 
Eq. (5) gives at once V=U, (7) 


which is here but a special case of the general proposition, easily proved by Eq. (1) and obvious 
physically, that if the potential is constant over part of the boundary of any region, and if the 


normal derivative vanishes over the remainder, the potential must be constant everywhere within 
the region. 


(2) Infinite half plane (Fig. 2) 


Choosing the origin at the division between I’, and I, we have for the general boundary conditions: 


v=f(x) :y=0, x>0; —0dv/dy= F(x), y=0, x<0. (8) 

We now transform to the w plane of Fig. (1) by means of: 
s=x+iy=w*=(§+in)?, (9) 
so that: r,:y=0, x>0>n=0; T,:y=0, x<0-§=0; f(x)f(#); (10) 


F(x)—2nF(— 1’). 


Putting these in Eq. (5), evaluating the integrals and returning to the z plane by means of Eq. (9), 
the result will be the solution to the original problem of Eq. (8). Furthermore, on inserting Eqs. (9) 
into the G of Eq. (4), the Green’s function for the half plane will be obtained in terms of (x, y). 

As a first example we shall take the case: 








=uUy:x>—-a@ 
filam; Pls) |. (11) 
=0 :x<-a@ 
Then: 
NVo 1 1 Uo f° #?+(no— 7)? 
v(&, n=— | + \aé—-— f no log dno 
wm Jo ((Eo—£)? +n? (Eo+&)? +77? nr Jo #+(no+n)? 
Uo 2+ (a+)? 2aé 
=et—| dant (+a? — 7?) log —————— —4né tan“! aaah (12) 
de +(a—n)? P+at—a? 
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On returning to the z plane by means of Eq. (9) the solution v will be obtained, explicitly in terms 
of (x, y). Thus for the negative x axis (¢=0) we find: 


a+ y—x\? 
[40 —x+(a’?+x) log (——) |: (13) 


Uo 
2r a-—-v¥-x 


v(x <0, y=0)=0+ 
As another example we may consider the case: 
f(x) =v0; F(x) = uo/x. (14) 


The solution for v then takes the form: 











no (” 1 1 uo f@ i + (no— 7)? 
v(t, 7) =— f | -+ Jag+~ — log dno 
wm Jo ((Eo—E)* +n? (E+E)? +7? r/o no 86 &+(mo+n)* 


=Vo9—Ugr+2uo tan E/n=vo9—U08, (15) 
where @ is the polar angle. 
This result can be readily checked by elementary means, as the system is simply the conjugate 
of that of a perfect radial potential distribution in which the equipotentials are circles about the 
origin and the streamlines or line of forces are the radii vectors. 


(3) Infinite quadrant with v or dv/dn given over parts of both axes 
This system (cf. Fig. 3) may be transformed to the normal quadrant distribution by the relation: 
w= (E+1n)?= (x+1y)? — x0", (16) 


so that the region where F is given goes into the +7 axis and the boundary where f is known goes 
into the +é axis. The resulting distribution will be given by: 








n ” 1 1 1 70 n 
v(&, n=" f Negetaey| + \ae.—— _ Fy (x0? — no”) }* log 
7 J0 (fo—E)* +7? (Eo+£)? +7? 2rdo (xo? — 0")! 


+ (no—n)? 1 - No #+(no—7n)? ” 
~ dy —Fe{ (not —xet)!} log dno _ (17) 
2+ (no+n)? 2a zo (o?—x0")! £+(notn)? 

where the meanings of F; and F2 are indicated in Fig. 3. 

(4) Infinite strip 


If dv/dy is given over the whole of y=), the Green’s function may be given directly in terms of 
(x, y) as: 





— {cosh 2(x—xo) /2b—cos r(y—yo)/2b} {cosh m(x—xo)/2b+cos r(y+yo)/2b} 18) 
=} log " 
* * feosh (x —x0)/2b—cos r(y+¥o)/2b} {cosh +(x —xo)/2b+cos +(y— yo) /25} 





When, however, dv/dy is given over only part of y=), the y axis may be set so as to go through 
the division point (cf. Fig. 4) and then the system can be mapped on the quadrant of Fig. 1 by the 


transformation : 
(E+1n)?=e**/>+1, (19) 
the region of x <0, y=b, being mapped on the segment 7=0, 0< <1. Furthermore, for §=0: 
dv/d§= (—2bn/a(n?+1))(d0/dy). (20) 


As an example we shall take the case: 


v=13 y=0; v=ve:y=b, x<0; dv/dy=0, y=b, x>0. (21) 
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Eq. (5) gives: 








V1 - 1 1 V2 ; 1 1 
v(é, )=— f + batt | me ~ d 
mt U(o-EP +n? (fot €)? +0? m Jo \(fo— EP +m? (Sot)? +7? 








Ve—V1 2n 
=y7,+ tan — 
T f+f—1 
=01+[(v2—v1)/m ] cos“!{ | ((e~**7/"+ cos ry/b)*?+sin® ry/b)'| —e-*2/>}, (22) 


(5) Semi-infinite strip (cf. Fig. 5) 


If v is given only along y=0 for x2 x» the transformation to the w plane of Fig. 1 is made by: 
(€+%in)?=cosh r(x+7y)/b—cosh rx9/b. (23) 

If v is given for all x>0, y=0, and x=0, OS y< yo, then: 
(€+in)?=cosh r(x+7y)/b—cos ryo/b. (24) 


Finally, if v is preassigned over the whole boundary except for y=b, x >xo, the transformation 
equation is: 
(€+in)?=cosh r(x+iy)/b+cosh rxo/b. (25) 


The rest of the procedure is exactly similar to that given for the other cases. 

If xo or yo in Eqs. (23) or (24) is zero the Green’s function may be obtained directly in terms of 
(x, y) by adding to the G of Eq. (18) a term exactly the same as given there except for a change in 
the sign in x9; and if xo in Eq. (25) is zero the above term should be subtracted from that in Eq. 
(18) to get G directly in terms of (x, y). 

As the problem of the closed rectangle would require elliptic integral transformations except when 
the point of separation between I; and I: coincides with one of the corners, and as the resulting 
integrals as in Eq. (5) would be in general intractable except by graphical treatment, this case will 
not be analyzed further. 


(6) Circle (cf. Fig. 6) 


The general boundary conditions may be written as: 
v=f(0) :0<0< Oo : dv/dr= F(8); O05 0< 2r. (26) 
The transformation to the infinite quadrant is accomplished by setting: 
(E+4n)?=1(20+2)/(zo—2) —ctn 49/2, (27) 


which, in polar coordinates, is equivalent to: 











2rro sin (0o— 8) re—r 
?—7= —ctn 69/2; 2§n= : (27a) 
re +r? —2rry cos (A9— 8) re+r—2rro cos (09—8) 
Ov —Aron Ov 
so that: (—) =— (—) ; (27b) 
0g t=0 nt— 27? ctn 60/2 +csc? 6/2 or r=r9 
As an example we shall take the case: f(@)=v; F(0) =u. (28) 





We then get: 
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2roun 7 No #+(no— 7)? 
v(€, 0) =v0— f log dno 
T o no'—2no* ctn 06/2+ csc? 69/2 #+(no+n)? 
2roto r+P— ps? w+? — a 2né 2né 
=Vo— tan-! — —tan~'! —————--- + tan"! ——————-- tan“! =| (29) 
8? — a’ 2nB 2na n+a?— rt+R-—# 
where: no'—2no? ctn 0o/2+ esc? 09/2= (ne? + a?) (n?+ 8). (29a) 


B. E1irHer T; oR Ty Consists oF Two SEPARATE SEGMENTS 


It will suffice for the present purpose to exhibit the transformation for the treatment of the single 
case of the infinite half plane, as all the others presented in Section A may be readily reduced to 
this problem. 

Thus if, as indicated in Fig. 7, the region: 0< x <x is that in which v is given whereas dv/dy is 
preassigned over the remainder of the x axis we may again transform to the w plane quadrant of 


Fig. 1 by the transormation : 
(¢+in)?=2/(xo—3), (30) 


the region x >xp mapping onto the 7 axis with 1£{ ~ and the negative x axis mapping on the 
rest of the positive » axis. The procedure now may be carried through in exactly the same manner 
as for the cases of Section A. 

When both I and I. consist of two separate segments, the infinite quadrant prototype of the w 
plane must be replaced by a finite rectangle. The given boundary may be mapped on the rectangle 
by means of an elliptic function transformation and then, as before, one may apply the Green’s 
function for this simpler system, namely : 


G(E, ) 
_ te {cosh #(&—2na— &)/b—cos r(n— no) /b} {cosh r(E—2na+ ko) /b—cos r(n— 70) /5} 
"a ” {cosh r(&—2na— &o)/b—cos (n+ 0) /b} {cosh r(E—2na+ &o)/b—cos r(n+ 0) /b\ 
where it has been supposed that the segments I’; are mapped on to: 0< <a, 7=0, b, and the segments 
rz to :0K Kb, E=0, a. 

Insofar as the boundary of any given geometrical region can be mapped on one of the simpler 
types treated above, and its parts I’; and I, do not contain more segments than permit a reasonably 
direct development of the Green’s function for the appropriate prototype boundary, the above 
method will give the solutions to the ‘“‘mixed’’ problems for those regions. Furthermore, it may be 
observed that with obvious changes Eq. (3) may also be applied to three dimensional systems, 
although its application to three dimensional problems suffers from the difficulty that there is no 
three dimensional analog of the conjugate function transformations, so that different Green’s 
functions have to be obtained explicitly for each different type of boundary. 

In addition to those problems occurring in electrical conduction in which the boundary values of 
the potential or its normal derivatives are usually constants, the mixed boundary value problem 
arises in the treatment of the irrotational motion of bodies in ideal liquids, in which the normal 
derivative is specified over the surface of the moving body and the potential itself at infinity. Heat 
conduction systems may also involve the mixed boundary value problem of the type treated here, 
if the surface radiation is arbitrarily specified and does not obey the Newtonian law of cooling. 
Finally, as suggested by Professor Bateman, it is likely that mixed boundary value problems in 
plane elastic theory, though governed by the biharmonic equation, may be reduced to equivalent 
problems in the logarithmic potential. The writer plans to return to this question at a later time. 

The writer wishes to thank Professor H. Bateman for reading the original manuscript of this 
paper and for his many helpful suggestions, and Dr. P. D. Foote of this Laboratory for the permission 
to publish it. 
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On Electric Charges of Disperse Systems Which Are the Seats of Physico-Chemical 
Reactions 


N. RAsHEvsky, Department of Physiology, University of Chicago 
(Received August 14, 1934) 


In a number of previous papers a mathematical study 
was made of small drops, which either absorb or produce 
some substances by virtue of physico-chemical reactions. 
This creates in the drop and in the surrounding medium 
gradients of concentrations of corresponding substances. 
In this note it is pointed out that when the reacting sub- 
stances are ionized, these gradients will result in the ap- 


pearance of electric charges and electric fields. An estima- 
tion shows that they are quite comparable with electric 
charges, produced by other factors in disperse systems. The 
fundamental property of the charges here considered is 
that they disappear as soon as the reactions cease. Their 
importance in biological systems is pointed out. 





HE question of electric charges of disperse 

systems has been a matter of much study. 
A variety of factors are involved in their pro- 
duction, such as the Donnan equilibria, the 
selective adsorption of ions, etc. All those factors 
are however of static nature. We wish to point 
out here another factor, of a dynamic nature, 
which must be present whenever physico- 
chemical reactions take place within a small 
liquid particle and when these reactions involve 
ionized substances. In a number of previous 
papers,' we have studied such systems and seen, 


that those reactions produce gradients of concen- ° 


trations both within and without the system. 
Those gradients are determined by the rate of 
reaction, the diffusion constants and the permea- 
bility of the drop-surface. If the substances 
involved are ions, then the gradients of concen- 
trations of ions of different nature will be differ- 
ent, and this will result in the appearance of 
free electric charges, though of course the sum 
total of all charges in the system will be zero. 
Consider the simplest case, that of two kinds 
of monovalent ions of opposite sign, present in 
equal amounts per unit volume in the sur- 
rounding medium at an infinite distance from 
the drop. And let the numbers of ions of both 
kinds, which disappear per unit time per cm* 
due to the reactions, be the same. For no charges 
can be created by the reactions themselves. The 
appearance of electric charges will produce an 
electric field, given for a case of spherical 
symmetry by d¢/dr, ¢ being the electrostatic 
potential. This field will affect the diffusion 
process of the ions. If c, denotes the number of 


1'N. Rashevsky, Physics 1, 143 (1931), 2, 303 (1932). 


positive ions per cm* and a, their mobility then 
instead of the ordinary diffusion equation we 
have, in polar coordinates, for a quasi-stationary 
state, as studied before: 


{@c,§ 2 dex! dg fdc..' 
D.'( nnn ilps —+2c.4/r) 
dr’ 
for the inside of the drop, and 


r dr dr \ dr 
d’c,° 2 dcx! deg (dc.* 
D( +) pasre"( +0.) =0, (2) 
dr? 


r dr dr\ dr 
for the outside, and corresponding equations for 
c_, the concentrations of negative ions. 
Further we have Poisson’s equation : 


(d?¢/dr*)+2(dg/dr)/r= —(4re/K)(c,—c_), (3) 





—q=0, (1) 








« being the electronic charge, K the dielectric 
constant of the medium (it will in general be 
different one inside and outside of the drop). 

At the boundary, for r=, we have: 


dc, dy dc..* 
D,—-= (|: _ —) (c4°—c+)= Di, 


(4) 
dr dr 


which, together with similar expressions for c_‘ 
and c_*, and with the requirement of continuity 
of g and dg/dr at r=ro, determines the problem 
completely. 

The treatment of this problem presents even 
in this simplest case very great difficulties 
because of the non-linearity of the differential 
equations. We shall here restrict ourselves only 
to a rough evaluation of the charges, which 
may thus be produced, by considering the case 


33 





34 N. 


of D,'=D,°=D_'=D_*=~@ and by neglecting 
the effects of the electric field on the diffusion 
process. We then have a uniform concentration 
of ions both inside and outside of the drop, and 
only a discontinuity at the surface. The concen- 
trations c, and c_ inside are now determined by: 
the condition that for a quasi-stationary state 
the amounts diffusing through the surface are 
equal to the amounts destroyed by the reactions. 


Hence: 
4arPhs(co—C4) =4 279°¢/3, 

and (5) 
4rrP?h_(co —_ c_) = 4rr,'q/3, 


since we assumed c,°=c_°=co, and q,=g_=q. 
Eqs. (5) give: 


C+ =Co— ro ‘Sh., c_=Co— ro, : = (6) 


Hence the total charge of the drop: 


4 4 1 1 
B= aardi(e.—)= carta ——— . (7) 

3 9 hk. he 
One point must be discussed in connection 
with this argument. We assumed the outside 
concentrations of both ions the same, and the 
rates of their absorption equal. Yet we find, 
after the quasi-stationary state is reached, an 
excess of one of the ions inside, without any 
change outside. The apparent paradoxon is due 
to our assumption of an infinite supply of both 
kinds of ions in the outside medium, so that 
drawing a finite amount from it does not change 
co. Actually there is an equal but opposite charge 

—E in the surrounding medium. 
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Unfortunately there are almost no data avail- 
able in the biological literature, which would 
give us / in absolute units. At any rate there is 
good reason to believe, that for most ions 
h<10-* cm: min.—'. For 4X10-* cm-min.“ is a 
plausible value for oxygen,? and the cell mem- 
brane is very highly permeable for the latter. 
Taking g~10-* g-cm-*- min.-'~10"* jons-cm=* 
-min.—', and assuming that /, differs only by 
1 percent from h_, we find for a cell with r»>~ 10-3, 
and h~10-> cm-min.—', a total charge of the 
order of 10’e. Under other conditions it may be 
much larger. A colloidal particle has usually a 
charge of about 10” ions per cm? of surface,* 
which for 7»~10-* cm would give for the total 
charge about 10’e also. 

We thus see, that the dynamic charges con- 
sidered here are at least of the same order of 
magnitude as the ordinary static charges of 
colloids. Their peculiar property is that they 
vanish with vanishing g, that is when the 
reactions for some reasons stop. 

This may have a particular bearing on charges 
of living cells, which are definitely connected 
with the life of the latter. If h,>h_ then the 
cell will be positive and a local increase of both 
h, and h_, that is a local increase of general 
permeability will produce a local negativity. 
More detailed considerations of the biological 
implications of the above will be given elsewhere. 


2 N. Rashevsky, Protoplasma 20, 125 (1933). 
3A. Gyemant, Grundzuge der Kolloidphysik, Vieweg, 
(1925). 
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The Mechanism of Division of Small Liquid Systems Which Are the Seats of 
Physico-Chemical Reactions 


II. Instability of Two-Phase Systems 


N. RaAsHEvsky, Department of Physiology, University of Chicago 
(Received October 19, 1934) 


In continuation of a previous paper, the same mathematical method is applied to a system, 
consisting of two spherical concentric drops, under the assumption, that a reaction takes place 
only in the outer larger drop, but not in the inner one. It is found, however, that even in that 
case the inner drop may become unstable for infinitesimal deformations, and spontaneously 


elongate. 





N a previous paper! we have investigated the 
mechanical stability of a drop under the 

action of such forces as result from the presence 
of concentration gradients, which gradients in 
their turn are produced by physico-chemical 
reactions, taking place inside the drop. We 
propose now to extend this investigation to more 
complex systems, consisting of two phases. In 
this paper we shall confine ourselves to the 
particular case of a system, consisting of two 
concentric spherical drops, of which only the 
external one is the seat of reactions, and we 
shall study the mechanical stability of the inner 
drop in such a system. On one hand this is one 
of the simplest cases of a two-phase system. On 
the other hand, as will be discussed elsewhere, 
this particular problem has some very interesting 
biological aspects. 

The mathematical method used will be the 
same as before. We consider first an undeformed 
system. Let R; and Rez be the radii of the inner 
and the outer drops, respectively. Let Do, D; 
and D, denote the coefficient of diffusion in the 
inner drop, the outer drop and the surrounding 
medium, correspondingly. Let h; be the per- 
meability of the surface of the inner drop, /»— 
that of the outer ; and let g be the rate of reaction 
per cm’ per minute, in the outer drop. Finally 





gR2 q(R2? —R;’) qgR? 


c=cot is eae a 
Iie 6D; 3D; 
gR2 q(R?—?r*) @Ri’ 

C= Co +— + ———_- - -—_ - 
le 6D; 3D ir 
q(R2*®—R;*) 1 

C-= Co + —_———_ -.. 

3D. r 


' Rashevsky, Physics 5, 374 (1934). 


let co be the concentration at infinite distance 
from the system. 

We shall confine ourselves to the action of 
osmotic forces only. Other forces discussed in 
our previous paper! and depending on concen- 
tration gradients, are easily taken into account 
in the same manner as before. 

Denoting by ci, ce and c; the concentrations 
in the inner and outer phase and in the sur- 
rounding medium correspondingly, for a de- 
formed drop we have to consider the sum of the 
osmotic potential RTc,/M and of the osmotic 
pressure RT(ce—c,)M at the surface of the inner 
drop. But this sum is equal to RT@/M. Hence 
we need concern ourselves only with values of cs 
at the surface of the inner phase. 

For the case of spherical symmetry we have 
the differential equations: 

V7°=0; DN*c;:—q=0; V%c.=0, (1) 
with the boundary conditions: 
for r= R, : Dj(dc;/dr) = hy(c;—c°) = Do(de®/dr), ‘ 
for r= Re : Dj(dc;/dr) = he(c.—c;) = D.(de./dr), 
where c°®, c; and c, stand for concentrations in 
the inner and outer phase, and the surrounding 
medium correspondingly. 

This set of equations is satisfied by: 





gRy*(D;— Reh2) bs q( Re’ — Ry’) 





3D; Rihe 3D.Re 
gRi(Di—Rehe) (Ro —R:*) 
* (3) 
3DiRihe 3D.Ro 
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We consider a deformation of the inner drop 
given by 
R=R,+A; A=Apcos 28. (4) 


At the deformed boundary we have, neglecting 
all higher powers of A: 


Ceijei+a=Leije,; (Cdei/drjrss=—gA Dey 
[ei—c Jay4a=0; [de°/dr Jri+s=0. 
Hence: 
[Do(de°/dr) —hi(c;—c°) Jrysa=9, 


(6) 
[D(de;/dr) —hy(e;—c°) Jryza= —gA. 


If therefore we represent again,' the solution 
for the deformed drop by 


C1=C°+03*; Co=Cit+c2*; cs=c.+c3*, (7) 
and require that 
V*c,=0; DN*c2—q=0; V2e3= 0, 


and that moreover we have 


D)(de,/dr) = hy(ce—¢}) = D(dc2 dr) 
D;(dceo/dr) = he(c3—c2) = D,(de3/dr) 


at the boundary of the inner and the outer drop, 





ao 


> A.P,.(0)=0; 


0 


with 


> B,P,.(0)=qAo cos 26; 
0 
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correspondingly, we find that ¢,*, c* and c;* 
must satisfy the following differential equations 
and boundary conditions. 


V°c,*=0; V%co*=0; V2c3*=0. (8) 
At the boundary of the inner drop: 
Do(dey* /dr) —hy(c2* —cy*) =0, 0) 
Dj(dceo* /dr) —hy(c2* —cy*) =A. 
At the boundary of the outer drop: 
Dj(dceo*/dr) —he(c3* —co*) =0, (10) 


D(dc3*/dr) —he(c3* —_ c2*) =0. 
We use again dc/dr instead of dc/dv (v—normal), 
for they differ only by terms of the order of A’. 
Eqs. (8) are satisfied! by: 


co*= Danr"P,(8), 
0 


c*=S(bar"+par-"-)P,(6), (11) 
0 


c3*= > sar~""P,, (0). 
0 


Introducing this into (9) and (10) and re- 
membering (4), we find: 


@o 


A,= Dona,R,"" —hy(b,Ri"+p,R-"" —a,R,"), 


B,= Dinb,Ri"'—Di(n+1)p,Rr-" 
€,= Dinb,R2""'—D(n+1)p,R2-"" 
Ds. 


> C,P,.(—)=0; ¥D,P,(6)=0, (12) 
0 0 
——fhy(b,Ri"+p,R1-"'!—a,Ri"), 
(13) 


2—he(s,R2-""'—b,R2* — p,R2*"), 


= —D.(n+1)s,Ro-"* —ho(s, Ro" —b, Ro" — p,»R2""). 


The first, third and fourth of Eqs. (12) give 
A,= €,=D,=0, (14) 

while the second gives after a simple calculation : 
Bo= —qAo/3; Be=4qAo/3, (15) 


all other %, being equal to zero. 

Hence only do, bo, po, So, a2, be, pe and Se are 
different from zero. Since we are interested only 
in the variation of the concentration along the 





surface of the inner drop, we need to know only 
the values with the subscript 2, for the 0-th 
terms of (11) do not contain 6. We find four 
linear equations for the determination of de, be, 
pe and se by putting »=2 into (13) and intro- 
ducing that into (12). But since we are interested 
only in the value of c2., we need to know only 
be and pr. 

After somewhat elaborate, but quite ele- 
mentary calculations, we find: 

















MECHANISM 


b= —cE/(AE+BD); p2=CD/(AE+BD) (16) 


with 
A=2[(D;—Do)iRv+2D .D;:Ri]; 
B= (3D;+2D )hiR1r?+6D DiRr; 


"= —4gAo(2DoRithR?)/3; (17) 


D= (3D,.—2D;)heRo? —6D;D.R2*; 
E= 3[(D.- D)h2R2? —3D;D.Ro* }. 


The use of those exact expressions is very 
inconvenient. We observe, however, that with 
D;, D. and Dy being of the order of magnitude 
of 10-'—10-* cm?-min.~!; with h; and he of the 
order of 10-4 cm: min.~! g~10-* g-cm~*- min. 
and R; and R.~10-* cm?) the first terms of A, 
B, D and E are very small as compared with the 
second ones. While in C the reverse holds true. 
Neglecting the small terms and putting R.=akR,, 
we find upon substitution of (17) into (16) the 
very simple expressions : 


be=2gAo/3D;a°Ri1; po= —4qAoR1*/9D;. (18) 


The approximation is good however only when 
a’>1. 

Introducing this into the second Eq. (11) in 
which we put r=R,, we find for the coefficient 
of P2(@) the expression : 


(2qgR1A0/3D 0°) an (4gRiA0/9D;)). (19) 

For a®>1.5 this expression is negative for 
positive g, and positive for q<0. A positive 
coefficient of P2(@) means that for @=72/2, that 
is at the equator, c;* is smaller than for 6=0, 
that is at the poles. And since, according to 
(5), c; has the same value along the deformed 
surface, a positive sign of (19) means a larger 
value of c; at the poles than at the equator, in 
other words it means a stable equilibrium for a 
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spherical shape. For g>0 the sign of (19) is 
negative, and the drop is unstable under the 
action of osmotic forces. Taking into considera- 
tion also the capillary forces, which always act 
in the direction of stabilization and which vary 
like Ao/R,’,? we find! that for small R, the 
capillary forces prevail and the drop will be 
stable even for g>0; but for larger R, the 
osmotic forces, which vary according to (19) 
like AoRi, will prevail and the drop will become 
unstable and may eventually spontaneously 
divide. In a similar way like before we find for 
the critical value of the radius: 


R,*=[36yMD;a°/RTq(2a°—3)]!, (20) 


R being the gas constant, T the absolute temper- 
ature, M the molecular weight of the reacting 
substance and y the surface tension of the 
inner interface. 

We find here, so far as the osmotic forces are 
concerned, a close analogy to the previously 
studied simpler case. A spontaneous division of 
the inner phase of a two-phase system can be 
caused by osmotic forces created only by sub- 
stances which are produced in the system. The 
interesting and important point of this investi- 
gation is however, that reactions, going on in 
the external phase only, can produce a division 
of the inner phase, though the latter may be 
chemically inert and even impermeable to the 
reacting substances. For our formulae (18) and 
(19) do not contain either Dy or i; and hold 
therefore with the same approximation for 
Dy=h,=0. The biological importance and impli- 
cations of this result shall be discussed elsewhere. 

The approximation used in deriving (18) is 
no good as soon as R, approaches 10-? cm. For 
such large drops the more complex expression 
(17) must be used. 


2N. Rashevsky, Protoplasma 16, 387 (1932); 20, 125 
(1933); Physics 2, 303 (1932). 
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The Flexure Problem for Rectangular Beams with Slits 


Dwicut F. GunperR, University of Wisconsin and Colorado State College 
(Received October 1, 1934) 


Experimental and approximate theoretical studies, made 
at the United States Forest Products Laboratory, Madison, 
Wisconsin, show that the mean shearing stress in the neu- 
tral plane of a rectangular beam with symmetrical hori- 
zontal slits or checks extending along its lateral faces is less 
than that given by the usual engineering formulas; that is, 
there is a definite ‘‘two-beam” action which relieves the 
shearing stress in the portion of the beam weakened by the 
slits. This paper presents an exact solution of the flexure 
problem for a beam of such section and confirms the main 
features of the distribution of shearing stress found in the 


earlier investigation. The solution of the flexure problem 
requires the solution of Laplace’s equation within the plane 
region coincident with a cross section of the beam, subject 
to certain boundary conditions. This solution is obtained 
by using the Schwartz-Christoffel transformation of a 
polygon onto a half-plane and solving the transformed 
boundary value problem. The exact solution is set up as an 
infinite series. To simplify the calculations two approximate 
solutions are also obtained, one of which is used in treating 
a numerical example. 





INTRODUCTION 


N experimental and approximate theoretical 

study' made at the United States Forest 
Products Laboratory (Madison, Wisconsin) indi- 
cated that the mean vertical shearing stress was 
less across the weakened section of a rectangular 
beam with symmetrically placed slits at the 
horizontal central line than that stress given by 
the usual engineering formulas. The following 
paper gives a more complete mathematical 
derivation of this result, presenting a method 
for calculating the shearing stress anywhere in 
the section of such a beam. 


THE FLEXURE PROBLEM 


In the analysis of the elastic behavior of a 
beam, two types of stress must be determined. 
One is the tensile or compressive stress due to 
the bending moment at the section, the other is 
the shearing stress due to the relative sliding of 
adjacent fibers. For metal beams the second of 
these stresses is comparatively small and can 
ordinarily be neglected in discussing the failure 
of the beam, whereas, for wooden beams of 
certain sections such as the one to be considered 
here, this is not the case. 

As is customary, the problem will be solved 
for a uniform cantilever type of beam with 
cross section of the shape under consideration. 
This problem, called the flexure problem, was 


1 


J. A. Newlin, G. E. Heck and H. W. March, New 
Method of Calculating Longitudinal Shear in Wooden 
Beams. U. S. Forest Products Laboratory, Madison, 
Wisconsin. Engineering News Record 110, 594-596 (1933). 
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first. solved by Saint Venant? who gave the 
method for computing the tensile and com- 
pressive stresses for beams of any section and 
obtained also the equations for determining the 
shearing stresses. The determination of these 
shearing stresses, which requires the solution of 
Laplace’s equation V’¥=0 subject to the con- 


dition 
f 2d ad 3 C 
/ 2 x——H+Co 
— =< 


¥ 


=¢= (1) 


‘ 


on the boundary S of the section has not been 
accomplished except for a few special sections 
such as the rectangle, circle, et cetera.* Recently 
W. M. Shepherd‘ gave the solution for beams of 
circular section with one or more radial slits 
extending inward from the boundary of the 
section.'In (1) uw, o and E are the usual constants 
of the material, Cy is an arbitrary constant, and 
the coordinate system is indicated on Fig. 1.° 
In terms of the function y, the shearing stresses 
car. readily be shown to be, 

} @ 


W dy | 


EI ay EI 

2 Saint Venant, J. de Math. (Liouville), Ser. 2 1, 1856. 

3A. E. H. Love, Theory of Elasticity, pp. 335-338, 4th 
Ed., 1927. 

4Shepherd, Proc. Roy. Soc. 
(1932). 

5 The equation ¢o=(W/ul)[@—x*/6], where @ is conju- 
gate to y, defines the relation between the function y of 
this paper and the function ¢o of Love’s general theory. 
See also, Griffith and Taylor, Tech. Rep. British Advisory 
Committee for Aeronautics 3, 963-969 (1917-18). Reference 
3, page 331, Eq. (4). It should be noted that due to the 
symmetry of the section the twist 7 vanishes. 


dy y° 
———+pox* 
ox 2 


London A138, 607-634 





FLEXURE PROBLEM 
METHOD OF SOLUTION 


The method which will be used in solving this 
problem consists of mapping the section onto a 
half-plane and solving the related boundary 
value problem. 

The Schwartz®-Christoffel’? transformation, 


F tes 
s=Cif TW(e-ajerae+C, (3) 


transforms any polygon in the z-plane into the 
upper half of the ¢-plane, the boundary of the 
polygon transforming into the real axis of the 
plane. In (3) the a; are the points on the real 
axis of the ¢-plane into which the consecutive 
vertices of the polygon transform, the za; being 
the corresponding internal angles of the polygon. 
By determining the constants of (3) for the 
section at hand and resolving the transformation 
into its real and imaginary parts x and y, the 
boundary function ¢ can be transformed into its 
value along the £é-axis. 

It is further known that the solution of the 
two-dimensional Laplace equation over the upper 
half-plane subject to the boundary condition 
that ¥ is to be equal to ¢ along the £-axis is 
given by, 

psin ag(£)dé 


1 co 
Wo, a)=— ff , -@) 
1 J-x t?—2tp cos a+ p? 


where (p, a) are the polar, and (é, 7) the rec- 
tangular coordinates of a point in the half-plane. 
A restriction which is satisfied here, and one 
which is sufficient to assure the validity of (4) is 
that ¢() be representable by a Fourier integral. 
After transformation ¢ is readily seen to be a 
function of ~ only and its direct substitution 
into (4) determines y. Finally, using this value 
and the relations (4) and (2), the values of X, 
and Y, can be found immediately. 

The boundary function (1) is determined only 
by the fact that the stress at a point on the 
boundary must be in the direction of the bound- 
ary at that point. From the symmetry of the 
section, it is evident that the x-component of 
the shearing stress vanishes at the vertical 
central line, hence either half of the section may 





6 Schwartz, Gesammelte Werke, Vol. II, pp. 65-83. 
7 Christoffel, Annali di Matematica Pura ed Applicata, 
(2) 1, 95-103 (1867); 4, 1-9 (1871). 


FOR RECTANGULAR BEAMS 
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B' 





C' *C-l-a () alc ; 
Fic. 1. Beam section before and after transformation. 


be chosen for transformation; for definiteness 
let it be the right half. 

Since the transformation (3) allows the arbi- 
trary choice of the points corresponding to any 
three vertices, respect for proper order being 
maintained, let the points O, B, D (Fig. 1) be 
transformed into 0, 1, ”, respectively, on the 
t-axis of the ¢-plane. The transformed coordi- 
nates a and c, of A and C will be dependent on 
the dimensions of the beam, i.e., on the width 
2 m, the depth 2 n, and the thickness of the web 
2 d. The expressions for a and c are obtained by 
writing Eq. (8) in its inverted form in terms of 
theta functions and obtaining from the initial 
conditions two relations involving only a, c, m, 
n and d. Following this the two values a and c¢ 
are readily found from these relations which are 
a?Ry—a(ey+1)+Roc;=0, and c*(c;?7—1) =c,?—a’, 
where ¢,?(1—k) =1+k and Ro=sn [(m—d)/7°M, 
k’}. 

After determining C; by proper choice of the 
lower limit ((=0, z=d) and replacing C; by a 
more convenient but still arbitrary constant 
— 71, Eq. (3) becomes 


J cdg 
s—d=—rif . 
0 ((¢?—a*)(¢?—-1)(¢?—c)}} 





(5) 


where the singularities of the integrand are to be 
avoided in the usual manner by enclosing them 
in infinitesimal circles and choosing the upper 
halves of these circles as the path of integration. 

As is seen by the aid of Fig. 1, the boundary 
function ¢ is reduced to (Cy chosen as n?m/2) 
= —pox®/3E+n*m/2 along OAB, to ¢= —pox?/ 
3E+n’x/2 along BC and to ¢=0 along CD, 
similar values holding for OA’B’C’D. At this 
point it is desirable to consider the relative 
importance of the different terms of ¢. With the 
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largest value of wo/3E encountered, the term 
uox®/3E is never more than approximately five 
percent of the value of @ at any point of the 
boundary and is much less than this for all 
beams of common occurrence. It is thus seen 
that for practical purposes this term may be 
neglected, the boundary conditions reducing to 
g=n'’m/2 along OAB, to ¢=n*x/2 along BC 
and to ¢=0 along CD. 

After introduction of the values of ¢ obtained 
from (5) by resolving it into real and imaginary 
parts, Eq. (4) becomes, 


nm 
V(p, a)=— pQ sin adé 
2nr vo 
n2 c 
+— [ x08 sin adi, (6) 
Q2rvi 


where the integrals for symmetric intervals 
about £=0 have been combined and where 
Q=(& —2Ep cos a+p?)!+(#2+2Ep cos a+ p?)—. 
The substitution, 


ry=2—[(2—a®)(e—1)]} 
6=+4[(e-a)(2-1)}) 
(2—a?)!— (1)! 


L " (c?—a?)!+(c?—1)! 


+¥-$* 


t= ————, where ¢ 
k(¢?—6) 








reduces (5) to 





é‘ dt 
s—d= -euif : 
wo [(@—1)(k—1)]} 














GUNDER 


Comparison of particular values of z, ¢, and ¢ 
corresponding to the vertices of the section gives, 


22MK=n, 72MK’=m, (9) 


K and K’ being the complete elliptic integrals 
with moduli k and k’=(1—?)!, respectively. 
Equating real and imaginary parts after the 
transformation (8), it is seen that on BC 


t<1/k dt 
v= m— Mf , y=n. (10) 
1 L[(@—1)(1—k#)]} 


Also by proper manipulation, Eq. (6) can be 
reduced to, 


¥(p, a) 
n?74 f é tan! 2tp sin a/(p?—#) 
2r 4, [(#—a*)(#—1)(2—c) }! 


which is the formal solution of the boundary 
value problem at hand. In order to complete 
the solution of the flexure problem it is necessary 
to invert the transformation (5) and to thus 
compute dy/dx and dy/dy and hence, by direct 











dé, (11) 


‘substitution into (2), the shearing stresses X, 


and Y,. 

Due to the difficulties which arise in inte- 
grating (11) even after the integrand has been 
broken up into a series, two simpler methods are 
given here® to aid in obtaining practical results. 


First method of approximation 


Investigation shows that the constant c (see 
above), is always near unity in value. Further, 














[ 2) Foss ! the function ¢ is monotone decreasing in the 
[c?—a? }}+[c?-1]} interval 1< €<c and a straight line approxima- 
where + tion for @ in this interval is found to be very 
tome eee (8) satisfactory. Using this the integral is readily 
ki - evaluated and yields, 
n?’m 2psina sina (c?+ p*)?—4cp? cost? a ¢ | 2cp sin a 
¥(p, a)= tan" OL- tan"! 
22 p—1 2(c—1) (1-+p%)?—4p? cos?a c—1 e—e 
2psina) cosa p’ sin 2a p’ sin 2a 
—tan™ {tan —tan™ || (12) 
p’—1 c-—1 c?>— p* cos 2a 1— p* cos 2a 


Second method of approximation 


A second method which yields a simpler value for y consists in obtaining two solutions, one 
with ¢=n*m/2 in 1£ &<c and one with ¢=0 in this same interval, and using the mean of these 


’ For a more detailed discussion of these solutions see the writer's thesis (University of Wisconsin Library). 
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values for y. The above two values of ¢ are, respectively, its maximum and minimum values in the 
interval and the results are accurate for points in the neighborhood of the slits, the error in the 
resultant shearing stress being not more than five percent. The value of y obtained by this method is, 


Y= (Wot 2) /2, (13) 


nm 2p sin a n?m 2cp sin a 
where, yo= — tan" —, and y2.=——tan"! 
2a p>—1 2a p?—c 





INVERSION OF THE TRANSFORMATION 


By resolving into real and imaginary parts, Eq. (5) or its equivalent form (8), is to be solved for 
p and a@ in terms of x and y. This solution will be expressed parametrically. Eq. (8) can be rewritten 


after change of the lower limit, in the form . 


f dt “( “) weal ' (14) 
=——i y—— §-+—(e— 5), 
o (4-A)1-RA)P on 2) m , 


which after inversion becomes, 





y¥-° K LK’ 
{= —_——__= sal —(2y —n)+—(m—x), é|- sn(u+iv) (15) 
k(¢?—5) n m 


defining the new variables u and v. The left hand member is resolved into real and imaginary parts 
by ordinary algebraic methods while the right hand member is similarly treated by using the addition 
formula for the sm function. The equating of these parts gives in terms of p and a, 


— p+ (6+) p? cos 2a—ys sn(u, k)dn(v, k’) 
k{_6? —26p? cos 2a+ p* ] 1—sn?(v, k’)dn?(u, k) 

















(16) 
p?(6—y) sin 2a _ sn(v, k’)en(v, k’)cn(u, k)dn(u, k). 
=4 : 
k[_6? —26p? cos 2a+ p* ] 1—sn?(v, k’)dn?(u, k) 
Solution of these yields (left hand members called U and V) 

kV(6—y) kV6 

tan 2a= ; = (17) 
k?(U?+ V*)6+k(6+y7)U+y kV cos 2a+(1+kU) sin 2a 


which are the required expressions in parametric form. The solutions (17), when taken with any one 
of the values of y, (10), (12) or (13), give a complete solution of the flexure problem for this section. 

Since the values of the shearing stresses are in terms of the partial derivatives of y the following 
additional formulas are necessary for computation, the notation 0/dx, y indicating two formulas, 
one with x and one with y appearing throughout. 


OU AV K’'sn(u, k)sn(v, k’)cn(v, k’) 
cnneiliincasiy Ck’? —dn?(u, k){1+dn?(v, k’)}], 
dx dy m[1—sn*(v, k’)dn?(u, k) ? 


dU —dV_ K’'dn(v, k’)cn(u, k)dn(u, k) 
count - [cn?(v, k’) —k®sn?(u, k)sn?(v, k’) ], 
oy Ox m(1 —sn*(v, k’)dn?(u, k) 














k(6—y) OV aU 
| (wat 7) +40 649) 47-2867) — {2k°6U+k(6+~7)} V: 
0a 2 Ox, y Ox, y 


ax, y (G7)? V+ [R5(U2+ V2) +R(5-+7) Ut} 





’ 
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Ox, y 


aV 
, ba (1+RU) sin 2a———2V{(1+kU) cos 2a—kV sin 2a} 
p 











Oa aU 
—kV sin 2a 
Ox, y Ox, y 














Ox, y 
OWo n?mn [ 
So p(p?—1) COS @ 
Ox,y (p?—1)*+ 4p? sin? a Ox, y 
Ove n?mc x! [ Oa 





= p(p?—c*) cosa 
Ox,y (p?—c*)?+4c*p? sin? al 


ay If Ovo dpe 
= + 
Ox,y 2L0x,y Ox,y 


Ox, y 





2plkV cos 2a+(1+kU) sin 2a} 





’ 


Op 
—(p?+1) sin rage 
Ox, Vv 


—(p?+c?) sin a— 


Op 
Ox, y 


| and X, and Y, frem Eqs. (2). 


The indeterminate forms of the above formulas which arise at certain boundary points of the 


CONCLUSIONS 


The above formulas have been used to com- 
pute the shearing stresses in a test piece whose 
dimensions were, m=1.25, n=2.25, and d=0.25, 
whence the values, k=0.6216, k’ =0.7833, a 
= 0.9524, c=1.0026, y =0.9825, 6=1.0281, were 
determined. 

Fig. 2 shows the distribution of shearing stress 
and the direction of this stress in the vicinity of 
the slits as calculated by the foregoing formulas. 
Those values omitted in the upper half of the 
figure have not been calculated since they do 
not affect the following discussions. The average 
maximum possible error, including that resulting 
from the omission of yox*/3E from the boundary 
function is not more than two percent. It is to 
be noted that this is probably much larger than 
the actual error. 

A very interesting phenomenon, which was 
discussed on the basis of experimental results 
combined with an approximate theoretical treat- 
ment by Newlin, Heck and March,! is shown 
in Fig. 3. This shows the vertical components of 
the shearing stress averaged over the width of 
the beam as this mean stress varies from the 
central line through the slits to the top of the 
beam. The graph very clearly proves the fact 
that, contrary to general belief, the mean 
vertical shear over the entire width of the beam 
is less at the central line through the slits than 





section are readily evaluated by elementary methods. 





it is just above or just below this line. In other 
words, when the checks are deep, the beam 
tends to act, not as one beam but more as two 
beams. This result confirms the work of the 
paper just quoted which found a pronounced 
“two-beam”’ action. This action was even greater 
than that indicated by the present treatment due 
to the fact that wood is not perfectly elastic in 
shear but is somewhat plastic with the result 
that the shear distortion in the plane of the slits 
before failure is greater than that which would 
occur in a perfectly elastic material. Calculations 
show that this two-beam action does not appear 
until the thickness of the web is reduced to less 
than half the width of the beam. 
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Fic. 2. Stress distribution in quarter section of beam. 
Fic. 3.,Mean vertical shearing stresses showing two- 
beam action. 
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A Determination of the Refractive Index of Vitreous Silica and the Calibration of 
Silica Refraction Thermometers Between 18° and — 200°C 


J. B. Austin ano R. H. H. Pierce, Jr., Research Laboratory, United States Steel Corporation, Kearny, N. J. 
(Received October 18, 1934) 


The variation with temperature of the refractive index of fused silica for the helium line 


VOLUME 6 


5877.2A has been determined between 18° and —200°C and the results have been applied to 
the calibration of vitreous silica refraction thermometers at low temperatures. These data have 
been correlated with measurements at high temperatures. It is found that the refractive index 
between — 200°C and 1000°C can be represented by the equation m= 1.45657 +-0.001553 ¢/100 


+0.00219 e9-33#/100, 





INTRODUCTION 


METHOD for calibrating vitreous silica 

refraction thermometers at high tempera- 
tures by means of a calculation based on the 
temperature variation of their refractive index 
has already been described.! The present paper 
gives the results of a determination of the 
refractive index of fused silica at low tempera- 
tures and extends the range of the calibration 
equation to — 200°C. The method employed was 
in principle the same as that used at high 
temperatures, namely, observation of the passage 
of interference fringes across a polished plate, in 
this case thermometer 311, when the temperature 
was changed by a definite number of degrees. 
The apparatus is shown in Fig. 1. The fused 
silica plate R was placed in a cavity in a small 
copper block C in which there was imbedded the 
junction of a copper-constantan thermocouple. 
The block itself rested on the bottom of a Pyrex 
tube P the upper end of which was closed by an 
optically flat window W sealed on with stopcock 
grease. This tube had two side arms, one leading 
to a vacuum pump capable of reducing the 
pressure to 0.002 mm Hg or better, the other 
providing an entrance, through a vacuum wax 
seal, for the thermocouple. The use of a vacuum 
in P served the double purpose of preventing 
the condensation of moisture on the surfaces 
of the plate and window, and of permitting the 
direct measurement of the change in absolute 
index of refraction, that is the index of refraction 
of the plate relative to vacuum, thus eliminating 
the correction for the change in refractive index 
of air with temperature. It also reduced the rate 


' Austin, Physics 3, 240 (1932). 





of cooling to a convenient value—approximately 
3 degrees a minute with the refrigerants em- 
ployed, so that satisfactory measurements could 
be obtained by the simple procedure of im- 
mersing tube P in the refrigerating bath and 
following the change in temperature of the plate 
with the thermocouple. On the other hand, the 
vacuum made it necessary to use great care in 
measuring the temperature because of the con- 
duction of heat along the wires. For example, in 
some preliminary observations in which the silica 
plate rested on top of the junction of the thermo- 
couple at the bottom of the tube, the inflow of 
heat along the wires was so much greater than 
the loss of heat across the vacuum to the silica 
that with the tube immersed in liquid air the 
thermocouple readings were high by almost 
100°C. The sudden admission of air reduced the 
observed temperature to very near its true value 
almost instantaneously. Errors from this source 
appeared, however, to be eliminated by the use 
of the copper block and by wrapping the thermo- 
couple around the block before bringing it to 
the side arm, thus insuring that at least 7 cm of 
the thermocouple were at the same temperature 
as the junction. 

The e.m.f. of the thermocouple was measured 
with a potentiometer with a sensitivity of five 
microvolts. The couple itself was calibrated by 
the Bureau of Standards immediately after the 
completion of the measurements; and as an 
additional check on the overall accuracy of the 
temperature measurements with the thermo- 
couple carefully purified mercury was substituted 
for the silica plate in the copper block and the 
melting point was determined from a heating 
curve. A flat lasting four minutes was obtained 
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Fic. 1. Cross section of apparatus. 


at —39.93°C, which is in satisfactory agreement 
with the value —39.87°C given in International 
Critical Tables (Vol. I, p. 53). It also confirms 
to some extent our belief-that the thermocouple 
wire was reasonably homogeneous and indicates 
that the unavoidable differences in thermal 
gradient in the couple during calibration and 
during use did not appreciably influence the 
readings. 

In making a determination the silica plate, 
copper block and thermocouple were set in 
place, tube P was sealed and evacuated, and 
the Dewar flask, containing the refrigerant was 
raised around it. The passage of the fringes on 
the plate was then observed and the temper- 
ature corresponding to each was recorded. Three 
refrigerants were used: an ice bath, an ether-CO, 
mixture, and liquid air, and observations were 
made over the temperature range obtainable 
with each of them. Since the initial temperature, 
that of the room, was different for each determi- 
nation all readings were corrected to the reference 
temperature 18°C, an adjustment which is easily 
made since the number of fringes passing per 
degree at room temperature is accurately known 
for the silica plate employed. 
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The final values, taken from a smooth curve 
which is the average of the several sets of data, 
are represented within 0.1 fringe by the equation : 


AN =2.6—0.131t— 1.90 10-42, (1) 


where AN is the number of fringes passing 
between 18° and £#C. From this empirical 
equation a similar expression can be obtained by 
dividing (1) by the thickness of the plate, 
0.465 cm, thus 


AN =(L)(5.59—0.281¢—4.09 X 10-4), (2) 


which, as in the case of the equation previously 
given for high temperatures,' we believe holds 
for any well-annealed thermometer used with 
the wavelength 5877.2A. A more general equa- 
tion applicable to any wavelength is obtained 
by multiplying (2) by \,/2 or 2.9386 10-5, or 


AN = (2L)/2,(16.45 —0.828¢— 12.0 
@ 10-42) x 10-5. (3) 


It should be noted that as the temperature is 
lowered vitreous silica becomes less and less 
sensitive as an indicator of temperature; thus, 
for thermometer 311, the limit of accuracy of 
observation, 0.1 fringe, corresponds to 0.6°C at 
room temperature and to approximately 2°C at 
— 200°C. Such a lack of precision is a definite 
disadvantage in many cases and, just as with 
most thermocouples at low temperatures, places 
a real limit to the usefulness of this method. 
If high precision is not required, however, this 
method provides a convenient and simple means 
of following temperature changes at low temper- 
atures in any system or apparatus which can be 
provided with a suitable window; and it is 
possible that a plate of some kind of glass with 
a higher sensitivity at low temperatures might 
be more suitable than silica as a thermometer. 

It is also possible to obtain from Eq. (1) an 
expression for the change in refractive index of 
silica in this temperature range by means of the 
relation 


An=AN),/2L—atn, (3a) 


where An is the change in refractive index 
between 18° and #°C, X, is the wavelength of 
the light used in vacuum, 5877.2A, a@ is the 
coefficient of thermal expansion of vitreous silica 
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and m is the refractive index at 18°C.? Values 
for \,, L and AN having already been given it is 
necessary only to evaluate a and ™, which need 
not be done with great accuracy since the 
quantity (atm) is but a small correction term. 
A review of the data collected by Sosman* shows 
that a can be represented with sufficient accuracy 
over the range 18° to —200°C by the equation 


a=10-*(0.44+2.67 X10-%), (4) 


while m, can be taken as 1.459. Substituting 
these values in (3a) we obtain 


An=1075(16.45 — 0.886¢— 15.90 X10~*#*) (5) 


for the change in m between 18°C and — 200°C. 

The values represented by Eq. (5) can be used 
to calculate the index of refraction at any 
temperature within the range specified if the 
absolute index of vitreous silica at 18°C for the 
wavelength 5877.2A is known. No direct deter- 
minations of this quantity appear to have been 
made with the requisite accuracy the only 
datum in the literature, that of Rinne,‘ being 
given to but five significant figures; but a value 
can be derived without difficulty by interpolation 
in the table of most probable values for 18°C 
at various frequencies as selected by Sosman.° 
The absolute index so derived is 1.45896, which, 
used with Eq. (5), yields the values for 0°C and 
lower found in Table I. The data for higher 
temperatures given in the table are calculated 
from An as given in the paper already cited! 
and are included to give an idea of the variation 
of m over the range —200°C to 1000°C. The 
data in this table can be satisfactorily repre- 
sented by the equation 


n = 1.45657 +0.001553¢/100+ 0.002 19e-% 33 #/100, 


as shown by the comparison in the table between 
the observed and computed values. 

The data in the table give no indication of a 
minimum in the refractive index such as is 
found in the density and such as would be 
expected if the equation of Gladstone and Dale,® 


* In the discussion of the high temperature calibration 
(reference 1, p. 247) Eq. (3a) is erroneously printed as 
An=,/2L(AN+ain,). The above form is correct. 

%Sosman, The Properties of Silica—Chemical Catalog 
Co., New York, 1927. 

‘Rinne, Neues Jahrb. Mineral. Beil. 39, 388 (1914). 
® Reference 3, p. 591. 
® Gladstone and Dale, Phil. Trans. p. 337 (1863). 
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TABLE I. The influence of temperature on the refractive index 
of vitreous silica for 5877.2A relative to vacuum. 











Nobs. Nobs. 
—Neale. —Neale, 
eC Nobs.* Neale.t 108 °C Nobs.* Neale.* X 106 
—200 1.45765 1.45770 5 500 1.46475 1.46475 0 
—150 1.45782 550 1.46546 
— 100 1.45808 1.45807 +1 600 1.46619 1.46619 | 0O 
— 50 1.45840 650 1.46691 
0 1.45880 1.45876 4 700 1.46764 1.46766 —2 
50 1.45923 750 1.46840 
100 1.45971 1.45970 1 800 1.46914 1.46915 —1 
150 1.46023 850 1.46990 
200 1.46079 1.46081 —2 900 1.47070 = 1.47066 4 
250 =1.46137 950 1.47144 
300 =1.46201 1.46204 —3 
350 =1.46268 
400 1.46335 1.46337 —2 


450 1.46404 








*based on ms = 1.45896 
+ from equation n = 1.45657 +0.001553 t/ 100 +-0.00219 ¢-0-33¢/100 


namely, (n—1)/d=constant=(n—1)v, where d 
is the density and v is the specific volume, were 
strictly valid in this temperature range. The 
temperature of minimum density seems to vary 
somewhat with the sample, or at least there is a 
considerable discrepancy in the temperature 
reported by the several investigators who have 
observed it, so that it is difficult to predict the 
magnitude of the effect upon —1 between the 
possible minimum and our lowest temperature, 
which was approximately —200°C. An average 
of the existing data for the density gives, 
however, a decrease between — 150° and — 200°C 
of close to 5 parts in 20,000. Since the constant 
in the Gladstone and Dale equation is about 
0.208, the corresponding change in should be 
about 5 parts per 100,000, which is equivalent 
to the reversal of a whole fringe and should have 
been easily observable. Inasmuch as no reversal 
was observed, the minimum, if it does exist at 
all, is not a pronounced one and the increase in 
n between the temperature of the possible 
minimum in our silica and — 200°C is less than 
5 parts per million; or it is possible that the 
temperature of the minimum may be close to 
—200°C. Values of the Gladstone and Dale 
constant for the wavelength 5877.2A have been 
calculated for the whole temperature range with 
results given in Table II. The constant shows a 
variation of only 3 percent over a temperature 
range of 1100°C indicating that the equation 
holds to a reasonable degree of accuracy for the 
change of refractive index of vitreous silica with 
temperature. 
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TABLE II. Temperature variation of the Gladstone and Dale 
constant (k=(n—1)v) for vitreous silica. 








eC s—i* 


— 200 0.45765 
— 100 -45808 
0 45880 
100 45971 
200 .46079 
300 46201 
400 46335 
500 46475 
600 .46619 
700 46764 
800 46914 
900 -47070 


v105* 





454.0 
453.9 
454.0 
454.0 
454.1 
454.1 
454.2 
454.3 
454.4 
454.4 
454.5 
454.6 








* taken from Table I. 
* taken from Sosman, reference 3, p. 361. 
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Fic. 2. Curves showing the effect of stressing the silica 
plate. Stress due to contraction of copper plate began at 
point S. 
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If, instead of taking m3 equal to 1.45896 as in 
Table I, we adopt Rinne’s determination, m3 
= 1.4592, a check on the consistency of his 
measurements at —64° and —160°C can be 
easily obtained. Thus, it is found that rounded 
values of m derived from Eq. (5) lead to n_¢s 
= 1.4586 and m_}69=1.4581, which are in exact 
agreement with Rinne’s data, showing that the 
change of » with temperature was the same in 
both cases. 

It is interesting to note in passing that in the 
first measurements in which the copper block 
was used, the cavity contracted more than had 
been allowed for and consequently began to 
squeeze the silica plate. The effect on the 
interference fringes was immediately evident as 
is shown by Fig. 2 in which the solid line repre- 
sents the observations for this run and the dashed 
curve is for the unstressed silica. In order to 
remove any strain which might have resulted 
from this treatment the silica plate was annealed 
at high temperatures before being used again. 


SUMMARY 


The change in refractive index of vitreous 
silica has been determined in the range 18° to 
—200°C by means of an interference method. 
The results are applied to the calibration of 
vitreous silica refraction thermometers in this 
temperature region. The number of fringes AN 
passing the reference mark on such a thermom- 
eter when the temperature is changed from 18° 
to @&C is given by the equation AN=2(L)/ 
(16.45 — 0.828¢— 12.0 10-*#) K10-5 where L 
is the thickness of the plate in cm and \, is the 
wavelength in vacuum of the light used. No 
minimum appears in the index of refraction at 
low temperature. 
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Magnetic Susceptibility of Some Alloys of ‘““Gamma-Brass” Structure 


Cyrit STANLEY SMITH* 
(Received September 26, 1934) 


The change of diamagnetic susceptibility with composition of a number of phases of gamma- 
brass structure (Cu-Zn, Cu-Cd, Ag-Zn, Ag—Cd) has been determined. In each case when the 
phase makes its first appearance with low electron ratio, it is not strongly diamagnetic, but the 
increase in diamagnetism is approximately lineally dependent on the number of free electrons, 
and the strongest diamagnetism occurs at the limit of the gamma-phase with the highest 


electron ratio possible for the structure. 





HE recent theoretical treatment by H. 

Jones! renders particularly interesting phys- 

ical measurements on that peculiar series of 

alloys of complex cubic structure known, from 

the first phase of the type examined, as ‘‘gamma- 
brass”’ structures.** 

Phases of gamma-brass structure occur fre- 
quently in alloys, but only at or around compo- 
sitions corresponding to a ratio of 21 valence 
electrons to 13 atoms, i.e., at compositions 
approximating Ag;Zns, CugAl,, Cuz:Sns, Fe;Zne, 
etc., the transition elements being considered as 
zero valent. This was pointed out by Westgren 
and Phragmen? and Westgren and Ekman* who 
extended the Hume-Rothery rule‘ to the gamma- 
alloys. The structure of all the phases is cubic. 
The unit cube contains 52 atoms and results from 
the removal of the center and corner atoms from 
a cubical group of 27 body-centered cubes and 
slight distortion of the remaining ones. The exact 
distribution of the different atoms in the different 
sets of equivalent positions varies with the 
particular alloy considered. The principal crystal- 
lographic determinations have been by Bradley 


*Research Laboratory, The American Brass Co., 
Waterbury, Conn., U. S. A. Work performed in The 
University Chemical Laboratories, Cambridge, England. 

** These are often referred to as “‘gamma structures,” a 
name that is apt to be misleading since many phases that 
are called gamma in the metallurgist’s nomenclature, are of 
totally different structures. This led Jones! to mistake the 
gamma-phase in the copper-cadmium system for the phase 
of gamma-brass structure, which in this system is the 
delta-phase. 

'H. Jones, Proc. Roy. Soc. A144, 225 (1934). 

2A. Westgren and G. Phragmen, Trans. Faraday Soc. 
25, 379 (1929). 


7A. Westgren and W. Ekman, K. Vetenskapakad. 


Arkiv. Kemi. etc. B10, No. 11 (1930). 
W. Ekman, Zeitz f. physik. Chemie B12, 57 (1931). 
‘W. Hume-Rothery, J. Inst. Metals 35, 295 (1926). 
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and Thewlis,°® 
Bernal.’ 

The results of Endo*® on a large series of 
carefully heat treated alloys, including, however, 
only two systems with a phase of gamma-brass 
structure, permitted Bernal® to make the gener- 
alization that the strong diamagnetism is a 
structural property and depends on the presence 
in the large unit cell of electron orbits including 
several atoms. 

The earlier measurements of Greulich reported 
by Weber!® showed a marked maximum in the 
diamagnetic susceptibility of cast copper-zinc 
alloys at about 32.8 percent copper, which was 
attributed to the compound CuZnp. 

The measurements of Spencer and John" and 
of Clifford,''* although they include systems in 
which phases of gamma-brass structure are 
present, are of little value in the present con- 
nection. As has so often been emphasized, it is 
useless to perform physical measurements on 
alloys unless they have been heat treated to reach 
equilibrium. A cast alloy or one that has been 
annealed for too short a time represents a quite 
indeterminate mixture of phases of irregular 
and variable composition, but suitable heat 
treatment will convert this to a structure in 
both chemical and physical equilibrium (except 
for grain size) and it is only on such material 


Bradley and Gregory,® and 


5 A. J. Bradley and J. Thewlis, Proc. Roy. Soc. Al12, 678 
(1926). 
( ot . J. Bradley and C. H. Gregory, Phil. Mag. 12, 143 
1931). 
7J. D. Bernal, Nature 122, 54 (1928). 
8H. Endo, Sci. Rep. Tohoku Imp. Univ. 14, 479 (1925). 
®J. D. Bernal, Trans. Faraday Soc. 25, 367 (1929). 
( 320) H. Weber (and K. Greulich), Ann. d. Physik 62, 666 
1920). 
pe ° F. Spencer and M. E. John, Proc. Roy. Soc. A116, 61 
1927). 
11a Q. C. Clifford, Phys. Rev. 26, 424 (1908). 
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that physical measurements should be made. 
None of the earlier workers have been primarily 
interested in gamma-phase, and their measure- 
ments on alloys were not sufficiently close to 
permit the behavior within the phase boundaries 
to be definitely established. 

The present investigation, prompted by Ber- 
nal’s discussion of the subject, was aimed 
primarily at determining the manner in which 
the susceptibility varies as the composition 
changes across the field of homogeneity. The 
systems copper-zinc, copper-cadmium, silver- 
zinc and silver-cadmium were selected. Of these, 
measurements have previously been reported 
only on the copper-zinc system. 


PREPARATION OF ALLOYS 


The copper-zinc alloys for the present work 
were made from pure cathode copper and a 
commercial grade of 99.99+ percent zinc. The 
other alloys were made from special zine or 
cadmium that had been redistilled from material 
containing 0.003 percent iron. Analysis of two of 
the copper-zinc alloys showed the presence of 
0.0065 and 0.0054 percent iron (alloys 31 and 41, 
respectively—Table I). A silver-zinc alloy (No. 
50—Table IV) contained 0.0012 percent, a 
copper-cadmium alloy (No. 61—Table IID), 
0.0002 percent, and a silver-cadmium alloy (No. 
67—Table V) contained 0.0003 percent iron. 
The silver used was of specially high purity but 
was not analyzed. Alloys of the desired com- 
position were made by first melting under a 
cover of powdered graphite or a flux consisting 
of LiCl and KCl the required quantity of 
the pure metals with intermediate alloys con- 
taining about 60 percent zinc or cadmium. 
The resulting ingots were reduced to coarse 
powder and then remelted under oxygen-free 
dry hydrogen* in vertical silica tubes, 0.55 

* Shimizu? has shown that vacuum melted alloys differ 
in susceptibility from alloys melted in contact with air. 
Goetz and Focke" did not find this effect in single crystals 
of bismuth alloys. It is obviously impossible to melt alloys 
containing zinc or cadmium in vacuum, and hydrogen was 
selected as the least harmful of the available gases to 
prevent too great volatilization. The volatility of cadmium 
and zinc result in a blanket of metal vapor which prevents 
much absorption of hydrogen,—indeed it is impossible to 
get blowholes due to hydrogen in brass, although hydrogen 
is a frequent cause of unsoundness in commercial castings 
of deoxidized copper alloys. 

"Y. Shimizu, Sci. Rep. Tohoku Imp. Univ. 21, 826 
(1932). 

18 A, Goetz and A. B. Focke, Phys. Rev. 45, 170 (1934). 
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cm bore, with the bottom end sealed. With an 
oxyhydrogen flame the silica tube was then 
melted off and sealed about 2 cm above the top 
of the specimen, the tube still being filled with 
hydrogen at atmospheric pressure. The speci- 
mens had an average grain diameter of about 3 
to 5 mm. 

The cast alloys, sealed in their silica tubes, 
were annealed for a long time to produce 
equilibrium structures and then were either 
quenched, air-cooled or very slowly cooled in 
the furnace. The last method was first used 
with the copper-zinc alloys and resulted in iron 
precipitation, so the other series were cooled 
more rapidly. Quenching was not possible with 
the copper-cadmium alloys on account of the 
extreme fragility of the alloys, but the silver 
alloys, which have the greatest mechanical 
strength of the series, could be quenched. The 
tubes containing the heat treated alloys were 
broken open with care and a sample for analysist 
chipped off from each end of the bar. This was 
done with a single sharp blow with a hammer 
and chisel which, after a little practice, gave a 
clean flat fracture without damaging the bar as 
a whole. . 


MEASUREMENT OF SUSCEPTIBILITY 


The magnetic volume susceptibility was deter- 
mined at 16 to 20°C by the Gouy method, by 
weighing small bars (0.554 cm) of the alloys 
with the lower end in a homogeneous magnetic 
field. The magnet, for the loan of which the 
author wishes to thank Dr. Cockroft of the 
Cavendish Laboratory, gave a maximum field 
strength of 18,500 gauss in a gap of about 
0.80.81 cm. The specimen was suspended 
by a fine nonmagnetic wire from a_ balance 
sensitive to 0.05 mg and protected from draughts 
by a tubular brass box with the necessary glass 
windows, the sides of the box being pierced with 
conical holes to fit closely around the poles of 

t The determination of copper was done by electrolysis, 
using a rotating cathode in a sulphuric-nitric acid solution 
of about 1 g of alloy. Silver was determined by titration 
with thiocyanate. Zinc or cadmium was taken by difference. 
In most cases there was a slight difference between the 
composition of top and bottom of the bar due to the 
distillation of volatile constituent, but this amounted to 
less than 0.2 percent in the copper-zinc system, 0.4 percent 
in the silver-cadmium system and 0.1 percent in the other 


systems. The reported composition is the average of the 
top and bottom samples. 
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the magnet. The box was extended to the base 
of the balance case, which was mounted on a 
sturdy shelf above the magnet and was furnished 
with means for fine adjustment for centering the 
specimen in the gap. The specimen was contained 
in a thin glass tube, and rested on two small 
projections made by pushing a point into the 
walls of the tube while hot. The top end of the 
tube had ears for attaching the suspension wire 
and at the bottom, 8 cm below the gap, was 
blown a small bulb which contained about 30 g 
mercury, the weight of which prevented the 
specimen from being thrown out of the field 
when strongly diamagnetic specimens were used. 
A correction for the specimen tube blank, 
amounting to less than 6 mg, was made to all 
readings. 

The weakest point in the method was the 
determination of field strength. This was done 


_by measuring the force acting on a column of 


water (the susceptibility of which was assumed 
to be —0.72X10~*) in a tube, with a cemented 
flat bottom, of internal diameter approximating 
that of the specimen. Since the specimens used 
were rather short (3 to 5 cm) the upper ends did 
not reach a point where the field strength was 
negligible and curves of (H?—H,?)! vs. magnet 
current were obtained for several specimen 
lengths. The true value of H was obtained by 
the use of a long column of water. The suscepti- 
bility calculations were made directly from the 
change of weight in the specimen and the change 
of weight of a column of water of the same length 
for the same current. Thus 


«-10°= —0-72Fd,?/Fyd?, 


where «x is the magnetic volume susceptibility ; 
F and F,, the forces acting on the specimen and 
on the column of water of the same length and 
of diameter d and d,, respectively, for the same 
current. F,, was read for each experiment from 
a large scale calibration curve. The current’ for 
the magnet (<15 amp., 110 volts) was obtained 
from a large battery of accumulators and was 
measured by determining the potential drop 
across a standard resistance with a Leeds and 
Northrup potentiometer. 

In all cases x was determined on each end of 
each specimen and at least five different field 
strengths varying from 10,000 to 18,000 gauss. 
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Fic. 1. Effect of heat treatment on magnetic susceptibility 
of some copper-zinc alloys. 


TABLE I. Magnetic susceptibility of some copper-zinc alloys 
in the cast condition. 








Composition (% Cu) 





Alloy No. Weight Atomic x(X 10°) 
26 38.32 38.99 — 2.84 
27 39.09 39.76 —2.70 
28 40.96 41.64 —2.72 
31 41.32 42.00 —1.79 
35 35.62 36.27 —5.03 
39 30.46 31.06 —5.27 
41 25.33 25.87 —3.44 








The resulting values of « plotted against 1/// 
gave straight lines and permitted x, to be 
determined, this figure being uninfluenced by 
ferromagnetic impurities. 

The density determinations were made by 
weighing in air and water, and the mass suscepti- 
bility, x, calculated in the usual way. 


EFFECT OF HEAT TREATMENT ON COPPER- 
Zinc ALLOoys 


In most cases the curves of x vs. 1/H were 
straight and almost horizontal, but the curves 
of the slowly cooled copper-zinc alloys which 
contained precipitated iron were steeply inclined. 
These showed the effect of heat treatment in an 
interesting way. The main series was annealed 
at 650°C and 300°C and very slowly cooled, a 
heat treatment intended to result in a state 
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closely approximating equilibrium at room 
temperature. This was satisfactory for the pure 
gamma-alloys, but those alloys containing the 
beta-phase were found to be paramagnetic and 
to have a susceptibility which varied enormously 
with field strength, an effect obviously due to 
the precipitation of free iron. The alloys not 
containing beta were diamagnetic and were 
uninfluenced by field strength, hence it is obvious 
that iron is almost entirely insoluble in beta at 
low temperatures but is soluble in the gamma- 
phase. This is to be expected since the iron-zinc 
system itself has a phase of gamma-brass 
structure and this should be at least partially 
isomorphous with other gamma-phases. The 
small amount of iron in beta can be rendered 
magnetically harmless by quenching from a 
temperature where it is in solution. Quenching 
from 600°C (and probably cooling at any rate 
other than an extremely slow one) holds the 
iron in solution. Fig. 1 shows clearly the effect 
of heat treatment on the susceptibility of two 
copper-zinc alloys, alloy 31 being in the beta 
+gamma-range and alloy 41 in the gamma- 
range. 

The results for the cast alloys (Table I) show 
no regular change with composition since they 
are not in equilibrium. Equally valueless are 
those measurements on the slowly cooled copper- 
zine alloys which contained iron precipitation 
which rendered the extrapolation to infinite 
field strength very uncertain. On account of 
greater purity and more rapid cooling, the 
alloys of the other systems studied were not 
paramagnetic and were practically uninfluenced 
by field strength. For all but the paramagnetic 
copper-zinc alloys the values of «,, are believed 
to be accurate to 0.01 10-*. 


RESULTS 


The results of the susceptibility measurements 
are given in Tables II to V and are plotted in 
Figs. 2 to 5 against electron ratio, the latter 
being the ratio of valence electrons to atoms, 
i.e., (Aa+Bb)/(A+B) where A and B are the 
atomic percentages* of elements of valency a 
and 6. The theoretical composition of the 

* For tables for the rapid interconversion of atomic and 
weight percentages in any binary system see Smith.'4 


¥C.S. Smith, Trans. Am. Inst. Min. Met. Eng. Con- 
tribution No. 60 (1933). 
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gamma-phases corresponds to an electron ratio 
1.62. This is not a critical point, but only the 
point where the structure is most symmetrical. 
It is immediately apparent from Figs. 2 to 5 
that the susceptibility of all the gamma-phases 
considered varies in the same way. Very little 
change in susceptibility occurs across the beta 
+gamma range, the susceptibility being of about 
the same order of magnitude as for the elements 
themselves. As the electron ratio increases across 
the range of the homogeneous gamma-phase the 
alloys become more strongly diamagnetic, and 
on the appearance of the next phase diamagnet- 
ism decreases once more, proportionately to the 


TABLE II. Magnetic susceptibility of gamma-phase in copper- 
zinc system. 








Composition Mag. Vol. Suscept. Den- 
(%Cu) —« (X 105) sity p, 
Quenched* g/cc, 


20°C 
8.19 


8.11 
8.08 


Mass Suscept. 


Alloy —x/p (X105) 
No. 


from An- 

600°C nealed Quenched 

(0.246) 

(0.296) 

(0.249) 

(0.159) 0.264 
0.249 0.316 
0.503 —_ 
0.629 —- 
0.775 0.767 
0.743 0.728 
0.715 0.725 
0.675 0.717 
0.604 0.589 
0.437 = 0.415 
0.420 — 


n- 
nealed* 


(2.01) 
(2.40) 
(2.01) 
(1.27) 
(2.01) 
4.04 
5.03 
6.16 
5.85 
5.64 
5.31 
4.70 
3.42 
3.37 


Weight 


46.65 
43.60 
41.30 
40.36 
39.39 
37.19 
35.44 
33.30 
32.29 
31.18 
30.37 
28.30 
25.26 
38.08 


Atomic 


47.36 
44.30 
41.98 
41.04 
40.06 
37.84 
36.17 
33.92 
32.90 
31.78 
30.97 
28.88 
25.80 
38.74 





0.209 
0.213 
0.228 


NN ee 
eww 


wrk UMmnUna 
nneue & © 


7 
7 
8 

+ 

2 
1 
7 
7 
6 
5 
2 








+ Alloys 29, 30, 36, 38 and 41 annealed 96 hrs. at 675°C, slowly cooled 
to and held 48 hr. at 500°C; slowly cooled to and held 120 hr. at 300°C; 
slowly cooled to room temperatures. 

Alloys 31 to 35, 37, 39 and 42 annealed 52 hr. at 675°C, slowly cooled 
to and held 96 hr. at 300°C; slowly cooled to room temperatures. 

Figures in ( ) are for alloys paramagnetic in weak fields due to free 
iron precipitation. The extrapolation to infinite field strength is rather 
uncertain. 2 

*a Alloys marked ® were annealed 141 hr. at 650°C and quenched in 
cold water; all other quenched alloys were annealed 5.2 hr. at 600°C 
before quenching. 


TABLE III. Magnetic susceptibility of delta-phase, copper- 
cadmium system (gamma-brass structure). 








Mass 
Suscept. 
—x/p (X 10°) 
Annealed 


0.247 
0.278 
0.296 
0.311 
0.378 
0.461 
0.592 
0.594 


Density 
p, #/cc, 
20°C 


9.04 
9.02 
9.04 
9.08 
8.87 
8.63 
8.76 
8.89 


Composition 

(YCu) . 
Weight Atomic 
34.25 47.95 
32.18 45.62 
30.18 43.33 
28.09 40.86 
26.16 38.52 
24.30 36.21 
22.22 33.56 
20.93 31.94 


—«( X 108) 
Annealedt 


2.23 
2.51 
2.67 
2.82 
3.35 
3.98 
5.19 
5.27 











*Annealed 268 hr. at 375°C, air-cooled. 
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TABLE IV. Magnetic susceptibility of gamma-phase, silver- 
sinc system. 














Mag. 
Vol. Mass 
Composition Suscept. Density Suscept. 
Alloy (%Ag) —«x (X10%) p, g/cc, —x/p (X10) 
No. Weight Atomic Annealed* 20°C Annealed 
47 56.93 44.48 3.14 8.87 0.354 
48 54.99 42.55 2.89 8.76 0.330 
49 53.93 41.50 2.84 8.75 0.325 


50 51.00 38.68 3.64 8.67 0.419 
43 50.94 38.62 3.54 8.65 0.409 
51 50.47 38.18 3.65 8.64 0.422 
52 49.11 36.90 4.02 8.61 0.467 
53 48.13 36.00 3.60 8.60 0.419 
54 46.74 34.72 2.88 8.53 0.338 








+ Annealed 65 hr. at 600°C, 96 hr. at 400°C, quenched. 
TABLE V. Magnetic susceptibility of gamma-phase, silver- 
cadmium system. 





Mag. 
Vol. Mass 
Composition Suscept. Density Suscept. 
Alloy (Sco Ag) —x (X10) p, g/cc, —K«/p (X10) 
No. Weight Atomic Annealed* 20°C Annealed 








64 35.20 36.15 3.32 9.60 0.346 


65 37.30 38.26 4.28 9.61 0.445 
66 35.43 36.38 3.53 9.64 0.367 
67 39.62 40.62 3.12 9.65 0.323 


68 40.67 41.66 2.88 9.73 0.296 


69 41.78 42.78 2.78 9.75 0.285 
70  ~=—-42.80 43.82 2.56 9.68 0.265 
71 43.22 44.24 238 9.73 0.245 

















* Annealed 239 hr. at 450°C, quenched in water at 18°C, 


amount of the latter. The diamagnetic suscepti- 
bility of the gamma-phase at the limiting 
composition corresponding to the lowest electron 
ratio is not much greater than that of the beta- 
phase (in the case of the silver-zinc system it is 
actually smaller), but there is a rapid increase 
with increasing number of electrons right up to 
the point where the epsilon phase appears. 
There are insufficient points to establish the 
type of curve exactly, but in each case they lie 
on a straight line within the accuracy of the 
measurements. These facts fit in perfectly with 
the theory of Jones.! 

The copper-zinc gamma-phase shows the 
strongest diamagnetism and the silver alloys 
the least. This is largely due to the fact that the 
gamma-structure remains stable to higher elec- 
tron ratios in the copper systems, and if it were 
possible to extend the limit of stability of the 
structure to still higher electron ratios, the dia- 
magnetism would undoubtedly become greater. 

The intersection of the different branches of 


SUSCEPTIBILITY 





OF ALLOYS 
the curves provide a means of determining with 
some accuracy the composition limits of the 
gamma-phases. The limits derived from Figs. 2 
to 5 compare favorably with those determined 
by metallographic or crystal structure determi- 
nations, although there are minor differences as 
can be seen in Table VI. A correction to the 
equilibrium diagrams cannot be made with 
certainty on the basis of the present results, 
however, since it is by no means certain that 
quenching retains the alloys in equilibrium, and 
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Fic. 2. The magnetic susceptibility of the gamma-phase in 
the copper-zinc system. 
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Fic. 3. Magnetic susceptibility of the delta-phase (gamma- 
brass structure) in the copper-cadmium system. 
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Fic. 4. Magnetic susceptibility of the gamma-phase in the 
silver-zinc system. 
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Fic. 5. Magnetic susceptibility of the gamma-phase in the 
silver-cadmium system. 


TABLE VI. ““Gamma”’-phase concentration limits. 








Atomic Percent Zinc or Cadmium 











System From Figs. 2—6 From Literature Authority 
CuZn 58.4-66.0 59.5-67.7 (380°) Owen and Pickup” 
57.8-67.1 (400°) Bauer and Hansen!* 
CuCd (6-phase) 60.1-67.2 57.4-63.9 (375°) Jenkins and Hanson"? 
AgZn 58.7-63.2 57.5—63.1 (380°) _ Owen and Pickup” 
AgCd 59.0-62.3 : 57.6-61.5 (450°) Durrant?! 
partial transformation might easily occur that ACKNOWLEDGMENTS 


would be undetected by any physical measure- 
ments. Susceptibility measurements actually 
made at high temperature, would provide an 
ideal method for the study of metallic equilibria 
and transformations. 

It was unfortunately not possible to determine 
the electrical conductivity of the alloys used for 
the magnetic measurements. The results of 
Puschin and Rjaschsky'® and Richards and 
Evans'* on the copper-zinc and copper-cadmium 
alloys, respectively, show, as would be expected, 
that both the conductivity and temperature 
coefficient increase with increasing electron ratio. 
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